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A brief critical summary of the principal contributions published on the theory 
of the subject is presented, followed by the report of a series of experiments designed 
to test the methods which have been proposed. It is concluded that while relative 
measurements can be made with sufficient accuracy, absolute values cannot be 
relied on as yet. The results emphasize the necessity for careful attention to such 
details as the construction of the camera, the preparation of the sample, and the 
characteristics of the x-ray film and of the microphotometer. 


I. INTRODUCTION 


HE x-ray method of determining particle-sizes depends on the fact that 

sub-microscopic crystals give rise to broad diffraction lines in the powder 
diagram. The apparatus used consists of a source of monochromatic x-rays, 
a slit-system to define the x-ray beam, and a cylindrical camera on the axis of 
which the sample to be studied is mounted. If the crystals composing the 
sample are small enough, the diffraction lines are broadened in excess of the 
value determined by the dimensions of the x-ray beam and of the sample. It 
is the purpose of this paper (1) to review the attempts which have been made 
to derive relationships between the particle-sizes and the breadths of the 
diffraction lines, since most of this work has been published abroad and may 
not be readily accessible, and (2) to give the results of our work on the ap- 
plication of these methods to particular cases. This work was undertaken not 
only to obtain answers to specific problems but also to find out as much as 
possible about the methods which have been proposed, their range of ap- 
plicability, their accuracy, and the most suitable techniques. The work is 
admittedly incomplete and several difficulties have remained unsolved, but 
it is believed that our experiences will be of value to others who are con- 
templating the use of these methods. 


II. SUMMARY OF THE THEORETICAL BACKGROUND 
In 1918, Scherrer! published the following equation relating line-breadth 
to particle-size : 
B = 2(In 2/7)'/*(A/A)(1/cos 0/2) + b (1) 


* Communication No. 103 from the Experimental Station of E. I. duPont de Nemours 
and Company. 

1 Scherrer, in Zsizmondy, Kolloidchemie 3rd Ed. pp. 387-409. 

2 The exact meaning of the term particle-size is given on page 59 in the explanation of 
von Laue’s formulas. 
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B is the line-breadth divided by R, the radius of the camera, A is the wave- 
length of the x-radiation, A is the particle-size, @ is the angle between the in- 
cident and the diffracted beams, and 3d is a constant depending on the di- 
mensions of the apparatus. The line-breadth is defined, in the derivation of 
all particle-size formulas, as: the ratio of the integral of the intensity across 
the diffraction line to the maximum intensity. This measurement would be 
difficult to make in practice and the quantity actually measured is the dis- 
tance between the two points on the line where the intensity is half the maxi- 
mum. The two definitions are, of course, equivalent only for diffraction lines 


of special shapes. 
A B Cc D 


Fig. 1. A, parallel radiation, transparent sample; B. parallel radiation, opaque sample; C, diver- 
gent radiation, transparent sample; D, divergent radiation, opaque sample. 
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Fig. 2. Line-breadth for infinitely large particles (after Brill). Abscissas: The half dif- 
fraction angle 6/2 in degrees. Ordinates: Relative line-breadth. A, divergent radiation, trans- 
parent sample; B, divergent radiation, opaque sample; C, parallel radiation, transparent sample. 


This equation was derived with the aid of the following assumptions: (1) 
that the crystal lattice of the sample is cubic, (2) that the particles are 
oriented at random, (3) that the absorption of x-rays in the sample can be 
neglected, (4) that the particles are similar in shape and equal in size, (5) 
that the lattice is undistorted, and (6) that the incident radiation is parallel. 
In practice the breadths, B, are plotted against 1/(cos 6/2); the slope of the 
resulting straight line is proportional to 1/A, while the intercept on the B 
axis is equal to b. 
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In addition to the limitation to crystals of the cubic system, the use of 
Scherrer’s formula has the serious drawback that parallel radiation and a 
totally transparent sample cannot be realized. Moreover, the best approxi- 
mations to these conditions which are experimentally feasible still leave much 
to be desired, as has been determined by Brill* and others. Scherrer expressed 
the breadth B of a diffraction line as the sum of two terms, the first a function 
of the angle @ and the size of the particles, the second a constant depending 
only on the dimensions of the sample and the slit system. For “infinitely” 
large particles the first term becomes zero and B is independent of 6. Now if 
the radiation is not parallel, the breadths‘ of the lines will decrease with 6; if 
the sample is not transparent the breadths will tend to increase with 0. See 
Figs. 1 and 2. Hence, these two defects will oppose each other and in favor- 
able cases the Scherrer condition may accidentally be realized. This appears 
to be the explanation of the excellent results obtained with the Scherrer 
equation in some cases, for example, Scherrer’s own determinations of the 
particle-size of colloidal gold.' 

In 1926, von Laue’ derived a formula which has been used as the basis for 
most of the subsequent work. His assumptions were the same as Scherrer's 
except that the lattice-structure of the crystals need not be cubic and that 
the radiation was assumed to be purely divergent. The expression is as fol- 


lows: 
wfr? ff] 
oo =) cos* — 
n\R 2 


wr cos? @\2 
1+ (< —) — 1 
nR 2 


which when solved for 7 becomes: 





n = (w/21)(B cos 0/2 — (1/B)(rr/ R)? cos 6/2) (2) 


where w=const. =0.55, r is the radius of the sample, R, the radius of the 
camera, and 7 is a quantity related to the particle-size by the following 


equation: 
N b,@\ 272 
n(hy, he, h3) = — | >(~ ) | (3) 
4r Mm; 


GO = Yohb/| SoAb;| and b; = [aay aac; 


Here \ is again the wave-length, h; are the Miller indices of the crystal 
planes concerned, 6; are the unit vectors of the reciprocal lattice, and m; are 
the numbers of unit lengths in the directions of the unit vectors a; of the 
elementary parallelopiped. The notation is that of von Laue except that 


where 


* Brill, Zeits. f. Krist. 68, 387 (1928). 

‘For purely divergent radiation and infinitely large particles the line-breadths 
B=(nr/R) cos 6/2, according to Eq. (2), page 59. 

* von Laue, Zeits. f. Krist. 64, 115 (1926). 
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he used x instead of 8. Moreover, his B had a different meaning.® For the cubic 
system Eq. (3) reduces to 


n = \/(4ram) = X/(47°A) (4) 


where A is again the particle-dimension. 

A condition which gives the limits within which the departure from purely 
divergent radiation can be neglected was also derived by von Laue as fol- 
lows: 


5/8w? cos 0/2(T/n)? = 1. (5) 


T is the Oeffnungswinkel, i.e., the pinhole diameter divided by the radius 
of the camera. 

So far all the theories had assumed that the sample consisted of particles 
of uniform size, obviously a fiction introduced for mathematical convenience. 
Patterson’ took over von Laue’s method and added the assumption that 
the sizes of the particles are distributed numerically according to the Max- 
well law, N,,)=y?e~?? where uw is a quantity proportional to the particle- 
size and N;,,) is the number of particles having the size corresponding to u. 
Patterson's calculation leads to the result that, if von Laue’s formula for the 
calculation of 7 is blindly applied to a material the particles of which are dis- 
tributed according to the Maxwell law, the value obtained will be approxi- 
mately equal to 4/7! where 7 is the mean of the values of » for the individual 
particles. The true mean particle-size cannot therefore be deduced without 
an exact knowledge of the distribution. Mark* objected to the assumption of 
a Maxwell law as favoring the large particles at the expense of the small ones 
and stated that a distribution law of the Gauss error type N,,) =e~“% 0? 
would seem a more reasonable choice. Without, giving the details of the cal- 
culation he added that if a material has the Gauss distribution, the result ob- 
tained by direct application of von Laue’s formula will be only slightly in 
error. ° 

Brill,’ and Brill and Pelzer®:° contributed three modifications of von 
Laue’s theory. These have the advantage that the assumptions more closely 
approximate actual experimental conditions. 

The first dealt with specimens which absorbed the x-rays appreciably. 
Assuming a solid cylindrical sample in which there was complete absorp- 
tion, Brill showed that for infinitely large particles the experimental line- 
breadths, B’, would be connected with the B used by von Laue by the equa- 
tion 


B = B’/sin? (6/2). (6) 


6 In von Laue’s Eq. (46), which has been quoted by several authors, the quantity B is ex- 
pressed in a different scale and must be multiplied by 2/(cos 6/2) before it can be applied to ex- 
perimental work. 

7 Patterson, Zeits. f. Krist. 66, 637 (1928). 

8 Mark, Trans. Far. Soc. 25. 387 (1929). 

® Brill and Pelzer, Zeits. f. Krist. 72, 398 (1929). 

10 Brill and Pelzer, Zeits. f. Krist. 74, 147 (1930). 
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His procedure was to correct the measured line-breadths by Eq. (6) and then 
to calculate particle-sizes by means of Eqs. (2) and (3). Finally since in any 
particular case absorption may not be complete in a thin layer at the surface, 
the calculated particle-sizes were plotted against sin? 0/2 and extrapolated to 
find the particle-size for sin? 6/2 =1 since at 6=180° absorption presumably 
has no effect. 

This treatment imposes two limitations upon the general application of 
the method. The accuracy of the final result obviously depends upon having 
one or more diffraction lines at large angles, which is often difficult to achieve. 
In the second place, it is not possible to get any information about the shape 
of the particles since the treatment tacitly assumes that all dimensions are 
equal. 

In addition to these limitations, it is perhaps worthwhile to point out 
some inconsistencies which may be sources of significant error. The correction 
for absorption embodied in Eq. (6) was developed for sharp lines, i.e., noth- 
ing was said. about the effect of absorption on the broadening due to particle- 
size in the derivation of this equation. Moreover, it was developed for parallel 
radiation while the von Laue theory was based on the use of purely divergent 
rays. The experimental work presented later will show that many of the 











Fig. 3. The hollow-sample method (after Brill). A, horizontal section through camera and 


sample (size of sample much exaggerated) ; B, double peak; C, single peak formed by very small 
particles. 


published results must be discounted for lack of attention to the definition 
of the incident x-ray beam. 

The second method® employed the sample in the form of a thin-walled 
hollow cylinder. The slit system had a single pinhole-opening on the circum- 
ference of the camera, which was placed as close as possible to the target of 
the tube so as to give divergent radiation. The rays encountered a greater 
thickness of material at the sides of the sample than in going through the 
tube and consequently the diffraction lines appeared double. See Fig. 3. The 
authors adapted von Laue’s theory to this arrangement and published a curve 
relating the separation of the peaks to the particle-size of the sample. When 
the particle-size is quite small the two peaks may fuse into one. The authors 
showed that the particle-size may be determined from a measurement of the 
line-breadth with the aid of a formula developed for this special case. The 
theory of this method is free from the inconsistencies of the first method; 
furthermore, it enables size-determinations to be made from a single diffrac- 
tion line, and hence the shape of the particles may be investigated. 
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The third method" is closely related to the second. The arrangement of 
the camera was the same, but the sample was made quite opaque. Materials 
whose absorption was not sufficient to accomplish this were coated thinly on 
a lead glass rod of suitable diameter. The theory postulates that only crystals 
in a thin layer of the material struck by the direct rays from the pinhole shall 
be effective in producing the diffraction lines. See Fig. 4. This means that in 
deriving the expression for the diffracted energy in a given line of the pattern, 
the integration is to be made over the arc BCD instead of completely around 
the circumference as in the method above. Unfortunately the authors in- 
tegrated over 180 degrees, i.e., the arc ACE, apparently to avoid a difficult 
mathematical problem. This amounts again to assuming parallel radiation 
in part of the theory. In Part III evidence will be presented to show that this 
inconsistency may be significant. The particle-size is to be calculated from a 
family of curves published by the authors. See Fig. 7. Values of b/r (the linear 
line-breadth divided by the radius of the sample) are plotted against @ and 
the corresponding values of 7 obtained by interpolation. This method, like 
the preceding one, enables a particle-size determination to be made from a 
single diffraction line and hence offers the possibility of investigating the 
shape of the crystals. In a later article Brill has studied this phase in more 
detail. 





Fig. 4. 


Ill. EXPERIMENTAL DETAILS 


The x-rays were produced by a Siegbahn-Hadding tube with inter- 
changeable targets excited at approximately 46 k.v. and carrying 13 to 15 
m.a. Copper and iron targets were used with their appropriate filters. Agfa 
x-ray film was used in all the work reported. Particular pains were taken to 
have the photographic conditions uniform. The films were developed for five 
minutes at 18°C using the elon-hydroquinone developer recommended for 
this emulsion; the trays were hand-rocked to ensure uniform development. 


Test of exposure-density relation and of microphotometer characteristics 


For the interpretation of the diffraction photographs it is almost essential 
that there be a linear relationship between the intensity of the x-rays pro- 
ducing a given deposit on the film and the deflection of the galvanometer 
when the spot is photometered, since any method of correcting for the want 
of linearity would be extremely tedious. To study this relationship, exposures 
were made on the type of film to be used in the subsequent work, the films 


" Brill, Zeits. f. Krist. 75, 217 (1930). 
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being held in a cassette 2 meters from the x-ray tube. On each film five or six 
exposures were made varied in suitable steps. Seven sets made when the tube 
was operating most satisfactorily were used to plot the curve shown in Fig. 
5. This procedure is only valid if time of exposure and intensity of radiation 
are strictly equivalent, i.e., y must be unity in the equation, Jf” =1, the so- 
called reciprocity relation. Bouwers” has shown that for x-rays y=1 to 
within one percent, so that the method is no doubt justified. 

Fig. 5 shows the relationship observed between relative exposures and 
deflections produced in the galvanometer. The curve closely approximates a 
straight line for low densities, and tests have shown that no appreciable error 
is introduced provided deflections greater than 35 mm, corresponding to 
densities“ greater than about 0.7, are avoided. The latter figure was obtained 
by photometering a calibrated neutral wedge since density determinations 
with the microphotometer alone were in error by as much as 50 percent due 
to the character of the illumination. 


100 [— 
DEFLECTIONS 





RELATIVE EXPOSURES. 





16) 100 200 


Fig. 5. Test of x-ray film and photometer characteristics. Abscissas: Relative exposures to x-rays. 
Ordinates: Deflections of galvanometer in mm. 


Tests of the Scherrer and von Laue methods with colloidal gold 


In order to get information as to the accuracy of these methods both as 
to relative and absolute values of particle-size and also as to the effect of the 
distribution of sizes in the sample, experiments were undertaken with col- 
loidal gold solutions of which the mean size and size-distribution could be 
determined by an independent method. 

A series of five sols was prepared following closely the method of Rinde.“ 
The nuclear sol was prepared by reducing chlorauric acid with phosphorus 
from water solution; sols of larger particle-size, by reducing with 33 percent 
hydrogen peroxide in the presence of a definite amount of the nuclear sol. 

” Bouwers, Zeits. f. Physik 14, 374 (1923). 

Density is defined as D = logo io/i where iy is the intensity of the incident light and 2, that 
of the transmitted. In photographic sensitometry the measurement is ordinarily made in diffuse 


light, which is not possible in the ordinary microphotometer. 
* Rinde, Dissertation, Upsala, 1928. In English. 
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By thus using each sol as a nuclear sol for the next larger, gold particles hav- 
ing any required size could be built up. The well-known precautions as to 
purity of ingredients and cleanliness of vessels were strictly adhered to. The 
sols were stabilized with gelatine, dialyzed with redistilled water for 24 hours 
or more to remove salts, and finally evaporated to dryness. Control experi- 
ments were performed to ensure that no step in this treatment altered the 
character of the sol. In each case a sample of the original sol was reserved for 
study in the ultra-centrifuge. 

The evaporated gold, which was somewhat plastic due to the gelatine, was 
formed into a cylindrical rod by packing it into a tiny glass tube and then 
extruding it about 3 mm from one end. The glass tube was then mounted on 
the axis of the camera. 

The camera used for these tests had a radius of 49.4 mm; the slit system 
was the conventional Debye-Scherrer arrangement with pinholes 1 mm in 
diameter. The diameters of the samples ranged between 0.4 and 0.7 mm. 











TABLE I. Particle-sizes of gold sols in angstroms. 








| C | C2 | C3 Cs C; 

Scherrer equation 47 62 | 81 

Laue equation 32 38.5 | 52 65 87 

Brill method Cu radiation} | 47 58 82 
Fe radiation) 30.5 42 | 

| 100 | 149 225 


Ultra-centrifuge 40 61 








Table I shows the results obtained with these sols. The first line gives 
sizes for three sols derived by means of the Scherrer equation, the values 
being found by the use of least squares. The second line gives results obtained 
by directly applying the von Laue equation which assumes no absorption in 
the sample. The next set of figures were found by the Brill method which 
allows for absorption in the sample but is otherwise identical with the pre- 
ceding. In the last line are figures derived from the published measurements 
on these sols by Rinde. These values were in good agreement with measure- 
ments made in this laboratory by Dr. J. B. Nichols.” Strictly speaking, the 
ultra-centrifuge determines the particle-size distribution on a weight basis, 
i.e., the distribution curve represents relative total weights of particles in 
the sample having a given diameter. The values given in Table I have been 
transformed to a number basis and represent the size of the most frequently 
occurring particles which is very close to the average size since the distribu- 
tion is nearly symmetrical. These values are of course not very accurate, but 
in view of the large discrepancy between them and the x-ray data this is 
unimportant. 

The source of this discrepancy is not obvious. The ultra-centrifuge method 
has received more intensive study than the x-ray, and has yielded particle- 
size determinations on gold sols and other materials in satisfactory agree- 
ment with those obtained by independent methods such as osmotic pressure 


For a brief description of the ultra-centrifuge and the method of particle-size determina- 
tion, see Nichols, Physics 1, 254 (1931). 
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and the ultra-microscope. At first sight, one would therefore be inclined to 
give greater weight to the values obtained by its use. On the other hand, the 
x-ray method is the only one which distinguishes between a single crystal 
and a crystalline aggregate. The fact that the disagreement in the results is 
relatively much worse in the sols of larger particle-size is perhaps evidence 
that there is some aggregation taking place. However, later work has shown 
how important is the geometrical definition of the x-ray beam, and conse- 
quently too great reliance cannot be placed on the x-ray results in Table I 
on account of the failure to provide either strictly parallel or purely diver- 
gent radiation. The agreement of the values obtained by the Scherrer form- 
ula with the other x-ray data must be partly fortuitous due to a favorable 
balance between absorption and divergent radiation as indicated on page 59. 

Additional conclusions which may be drawn from this work are, first, that, 
although the absolute values may be seriously in error, any of the three 
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Fig. 6. Size distribution of a colloidal gold solution. Abscissas: Particle-radius in angstroms. 
Ordinates: Relative number of particles. Note: The values for small particles are too high on 
account of a systematic error introduced by the microphotometer which has since been elim- 
inated. 


methods used is capable of arranging the series in the correct order of size 
with considerable precision. In the second place the particle-size distribution 
curves as obtained by Dr. Nichols for these sols, when placed on a number 
basis approximate quite closely the error type of curve and not at all, the 
Maxwell law. Fig. 6 shows this result for one sol. According to Mark, there- 
fore, the effect of the distribution of sizes on the value obtained by x-rays will 
be small. 

It might also be mentioned that this work furnishes direct proof—for the 
first time, it is believed—that a colloidal particle formed by the nuclear 
method does grow by extending the crystal lattice, as indeed has long been 
supposed. 


Tests of the Brill and Pelzer hollow-sample method 


The results reported in the preceding section indicated that insufficient 
care had been taken by the writer and by others in defining the character of 


16 In some preliminary experiments using samples of lower atomic weight, negative values 
of particle-size were obtained. The explanation for this is suggested by Fig. 2. 
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the x-ray beam. Accordingly a camera was built to give as nearly purely 
divergent radiation as possible. An examination of Eq. (5), above, shows that 
with cameras of reasonable dimensions only samples having extremely small 
particle-size will satisfy von Laue’s condition. For example, given a camera 
of 10 cm diameter with a pinhole 0.25 mm in diameter the largest particle- 
size which will satisfy Eq. (5) is about 60A. 

On account of the great difficulty of allowing mathematically for the 
finite size of the pinhole, the device of varying its size experimentally was 
resorted to. Three pinholes were constructed having diameters of 0.66, 0.46 
and 0.24 mm and arranged so as to lie accurately on the circumference of the 
camera. A sample of magnesium oxide was prepared by casting a thin film 
from a suspension in collodion, wrapping it around a copper wire and slip- 
ping the wire, which had been very lightly greased, out of the tube. Diffrac- 
tion photographs were then made with each of the pinholes without disturb- 
ing the sample. The films showed clear-cut double lines which gave excellent 
photometer records. The particle-size values, which are tabulated in Table 
II, are seen to vary considerably with the size of the pinhole. Brill and Pelzer® 
also used this method to determine the particle-diameter of magnesium oxide 
and reported the values given in the last column of the table. They did not 
give the dimensions of the apparatus used. 


TABLE II. Particle-size of magnesium oxide in angstroms. 








Brill & 


7” Pinhole 0 mm 
Indices aan com 0.46 mm 0.24 mm extrapolated Pelzer 
100 52 59 65 72 49 
110 57 61 81 88 57 
111 80 95 106 54 

















Attempts were made to apply the hollow-sample technique to titanium 
dioxide, colloidal silver and two different preparations of colloidal gold but 
without obtaining any satisfactory diffraction patterns. Titanium dioxide 
alone gave one line which showed a suggestion of the double peak. It was 
concluded that these materials are too absorbing even in such thin tubes and 
that the method can only be applied successfully to very transparent sub- 
stances. 


Tests of the “total absorption” method with colloidal platinum 


Four platinum black catalysts made by different procedures'’ were used 
to test this method. The samples were prepared by very lightly greasing a 
lead glass rod of appropriate diameter, dusting it with the colloidal metal 
and finally rolling the preparation on a glass plate. The photographs were 
made in the camera previously described, all three pinholes being used in 
most cases. Exposures of from 8 to 26 hours were required depending on the 
size of the pinhole. At least four diffraction lines were measurable in all 
cases. The results are displayed in Fig. 7. 


17 Taylor, Kistiakowsky and Perry, J. Phys. Chem. 34, 748 (1930). 
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In the case of Pt. Black No. 1, curves for the largest and smallest pinholes 
are plotted. There was again a strong dependence of the apparent particle- 
size on the diameter of the pinhole, as the curves and Table III show.The 
points fitted the theoretical curves well with the exception of the measure- 
ments for the first two diffraction lines, which, in every case, fell above the 
curves. The only reasonable explanation of this which can be suggested is 
that the layer giving rise to the diffraction lines at small angles was not suf- 
ficiently thin or opaque and hence gave lines of excessive breadth. 
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Fig. 7. Particle-size determination by the total absorption method. Dotted lines are theore- 
tical curves. Solid lines are experimental: Crosses—Pt. Black No. 4, Large Circles—Pt. Black 
No. 3, Squares—Pt. Black No. 2, Triangles—Pt. Black No. 1, largest pinhole, Small Circles 
Pt. Black No. 1, smallest pinhole. Abscissas: Diffraction angle @. 





TABLE III. Pt. Black No. 1. 








Pinhole diameter Particle-size 
0.66 mm 290A 
0.46 313 
0.24 430 
0 extrapolated over 500 
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Pt. Blacks Nos. 2 and 3 gave very diffuse diffraction lines indicating that 
they were much more finely divided than No. 1. When the data were plotted 
it was observed that they no longer fitted Brill’s curves, the lines being much 
too narrow at small angles. In order to vary the experimental conditions, a 
new platinum preparation was made and photographed with iron radiation; 
the previous results were all obtained with copper. The resulting curve, 
designated as Pt. Black No. 4, is consistent with those for Nos. 2 and 3. 

The results suggest that the theory may be in error because of the incon- 
sistency referred to on page 62 and in Fig. 4. Although the inaccuracy in 
placing the limits of integration would have little effect on the breadths of 
the lines at large angles, it might be quite important at small angles and 
would act in the direction required by the curves presented, i.e., the theore- 
tical breadths at small angles would be less. 

If the line-breadths at large angles are taken as more probably correct, 
the particle-sizes obtained with the smallest pinhole are, for Pt. Blacks Nos. 
2 and 3, about 60A; for Pt. Black No. 4, about 50A. 


Experiments with mixtures 


In order to obtain some idea experimentally of the effect of the distribu- 
tion of sizes it had been planned to make up samples having various artificial 
size-distributions which could be determined with the ultra-centrifuge and 
hence correlated with the x-ray findings. It was not possible to carry out 
this program; however two preliminary experiments will be reported for 
what they are worth. Two mixtures of Pt. Blacks Nos. 1 and 3 were prepared 
and particle-size determinations made by the total absorption method using 
the smallest pinhole. The results were as follows. 


Percent of Pt. Blk. No. 1 Apparent particle-size with 
in mixture smallest pinhole 
0% 60A 
20 186 
50 240 
100 430 


The diffraction lines due to the mixtures were abnormal in shape being 
obviously the result of the superposition of two lines of the shape ordinarily 
observed. This is evidence that each of the Pt. Blacks used is relatively uni- 
form in size. It also suggests that the shape of a diffraction line can be used 
as a qualitative criterion of the uniformity of an unknown sample. 


IV.. CONCLUDING REMARKS 


The experiments reported indicate that in many cases the results obtained 
may err by a factor of 2 in absolute value though it appears that relative 
estimates of particle-size on similar materials may have considerable signifi- 
cance even when some of the refinements are dispensed with. It is probable 
that the question of particle-size distribution can be overlooked at the 
present time although determinations on materials which have been pre- 
pared under controlled conditions can obviously be relied on to a greater de- 
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gree than when the history of the material is quite unknown. This paper has 
not touched on the problem of lattice distortion, although most of the mater- 
ials worked with have been prepared under such conditions that this source 
of trouble might be expected to be minimized. 

Further work, both theoretical and experimental, is obviously indicated, 
especially with regard to the total absorption method of Brill and Pelzer 
which is preferable to the others, particularly if the theoretical treatment can 
be carried through more rigorously. 

The writer wishes to express his appreciation of the many helpful dis- 
cussions with Dr. A. L. Patterson of the Johnson Foundation, University of 
Pennsylvania, and with Dr. A. W. Kenney of this laboratory who suggested 
the study. Thanks are also due Dr. J. B. Nichols for his ultra-centrifuge work 
and Mr. H. Aughey for skillful technical assistance. 
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Thickness Vibrations of Piezoelectric Oscillating Crystals 


By Issac KoGa 
Tokyo University of Engineering 
(Received April 4, 1932) 

For several years, it has been well known that, in the natural vibration of a 
piezoelectric oscillating crystal, there is a certain period proportional only to the 
thickness, but, so far as we are aware, there has been no concrete explanation for 
this vibration nor how its period is connected with the physical constants of the 
medium. The present paper shows that such a vibration, named thickness vibration 
for brevity, is due to the standing wave produced by interference of plane waves 
incident to and reflected from the plane boundary surfaces of the medium, and veri- 
fies the theoretical results by several examples. There are three normal modes in the 
thickness vibration, and the corresponding frequencies are given in first approximation 


9 


by a representative formula: P/2%=(q/2a) (c/p)''*, where g is any integer, a the 
thickness of the plate, p the density and c¢ a certain adiabatic elastic constant de- 
pending upon the orientation of the plate with respect to the crystallographic axes 
of the medium. It is not always possible to excite all normal modes of vibration piezo- 
electrically. The electrically measured natural frequency of the thickness vibration 
is always a little higher than that calculated by the formula. Adiabatic elastic con- 
stants of piezoelectric crystals may be determined in first approximation by measure- 
ment of the frequencies of thickness vibrations of plates prepared from given crystals. 


I. INTRODUCTION 


S IS well known, a thin piezoelectric crystal plate has in general many 

modes of vibration, but the shortest period of fundamental vibration is 
always dependent only upon the thickness of the plate. We shall name this 
vibration the thickenss vibration. 

Although many papers have been published upon the vibration of piezo- 
electric crystal, none has yet pointed out clearly, as far as we are aware, that 
the thickness vibration of thin crystal plate is a mere result of standing 
waves produced by the interference of plane waves, the wave surfaces of 
which are parallel to the boundary surfaces. 


II. EQUATIONS OF MoTIoNs IN CRYSTALLINE MEDIA 


The fundamental period of the thickness vibration of a very thin crystal 
plate is dependent almost entirely upon the thickness of the plate, so that 
we may consider the extent of the surfaces to be infinite. 

The free vibration of a thin crystal plate having two parallel plane bound- 
ary surfaces is due to the interference of plane waves incident to and re- 
flected from those boundary surfaces. In other words, in the thickness vibra- 
tion of the crystal plate bounded by the two infinite parallel planes x =0 and 
x =a, the displacement at any point in the medium is dependent only upon 
x and ¢, and not upon y and z. Therefore the general equations of motions in 
a crystalline medium, 
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OX ,/Ox + OX,/dy + O0X,/dz = pd*u/dt’, ) 
OY ./dx + OY,/dy + OY./dz = pd*r/dt? |} (1) 
OZ ,/Ax + AZ,/dy + IZ,/dz = pd?w/dt?, | 

where 


> 
& 
| 


= Cy1€rr + C2 yy + €13€zz + C14 yz + C15€zz + C16€ ry 


Xy = (1€r2 + C62€ yy + C63€22 + C64€ yz + Coszz + C66€ ryy 


Chk = Ckhy 
(2) 
Ou Ov Ow 
Crz =) Cyy =_-—) Czz le ) | 
Ox dy Oz 
Ow + Ov Ou Ow Ov Ou 
= = —? Ge O— pe) Cop ®@ —- +; 
Oy Oz Oz Ox Ox dy } 


u, v, w, displacement in the direction of x, y, z respectively, p, density of the 
medium, reduce to the very simple forms: 
OX, O7u oY, Oy OZ, Ow \ 


= g—> —_= ; —- = —y , (3 
ax “ae ax "ar ax "ar § o 


or further 
(0?/0x")(c\\u + C16v + C15W) 
(07/O.x) (Cig + Cog + CrgW) = pd*v/dt?, f (4) 


pda*u/dt?, —) 


(0?/Ax*) (C15 + C560 + C556) = pd? w/dl?. 
Now let /, m, n and c be so chosen that 
ley, + mei + ne15 = le 
lcig + meg + NCsg = Mc 7+ (5) 


leis of MC56 + NC55 = NC 
and put 
—& = lu + mv + nw, (6) 


then, multiplying the three equations of (4) by /, m, n and adding together, 
we have in the most compact form: 


0°E/dl? = (c/p)(0*E/0x"). (7) 
Eliminating /, m, n between the three equations of (5), we have 
cu ¢ C16 C15 
C16 Cogs — C C56 = 0) (8) 
| Cis C56 fu = € I 


a cubic equation determining c; let its roots, which are proved to be neces- 
sarily real and positive,’ be c, C2, cz. When c; is substituted in any two of the 
foregoing equations, the ratios of 1:m:n can be derived; let them be denoted 


1 Lord Kelvin, Baltimore Lectures, London, 1904. 
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/ 


by /;:m,:m, and suppose the corresponding value of & to be & with similar 
expressions for the other values of c. Then &, £2, and £3 are always proved to 
be perpendicular to each others. 

The solution of (7) is a very simple matter. To determine the normal 
mode of vibration we must assume that & varies as a harmonic function of 
the time ¢ 

Eo e'Pt, (9) 


Then as a function of x, — must satisfy 


d*§/dx* + (p*p/c)ét = 0 (10) 
of which the complete integral is 
& = A cos [px(p/c)'/?] + Bsin | px(p/c)'/?| (11) 


where A and B are independent of x. 
Now when both boundary surfaces are free from tractions, 


X,=Y,=2Z,=0 or 0&/dx = 0 at x = 0 and x = 4a, (12) 
we get 
B = 0 and sin [pa(p/c)'/?] = 0 (13) 
from which 
p/2m = (q/2a)(c/p)'!? (14) 


q being integral. 
Accordingly, the normal mode is given by equation of the form 


£ = A cos (qrx/a)-e'”" (15) 


where A is an arbitrary constant, which may be determined in the usual 
manner, when the initial value of & is given. 

Similarly if the medium be bounded by two parallel planes y=0 and y =a, 
we may replace (5), (8) and (15) by 


lege + moog + Nye = Ic; leap + Moe + NCoy = me; Ie4g + Moy + NCgy = Ne. (16) 


| C66 — € C26 C46 
| C26 Cea = € Co4 = () (17) 
| 
C46 C24 Cau — CI 
— = A cos (qry/a)-e'?! (18) 


and if the medium be bounded by z=0 and z=a, by 


less + meg5 + C35 = Ie; e435 + meg + C34 = me; 1035 + Me34+ NC33 = Ne. (19) 


Cs5 — € C45 C35 | 
C45 Ca — € C33 0~= ||| = O (20) 
C35 C34 C33 — € | 
— = A cos (qrz/a)-e'?'. (21) 
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It is well known that we can excite one or more of the free thickness 
vibrations by applying an alternating electric field of suitable frequency in 
the direction normal to the boundary surfaces of crystal, provided the crys- 
tal plate is prepared so that it is piezoelectrically active; and also that it is 
sometimes possible to sustain a continuous oscillation by properly employing 
the thermionic tube circuit. The ease of maintenance depends considerable 
upon the medium and the orientation of the boundary surfaces. 





III. THICKNESS VIBRATION OF QUARTZ PLATE 


Of the piezoelectric crystals, quartz is the most familiar to high-frequency 
engineers, so that we may first verify the theory with this medium. In 
quartz, if the coordinate axes be chosen as in Fig. 1, the adiabatic elastic 
constants? and density are as follows: 


Cu = 85.45 X 10!’ dynes/cm?, C2 = 7.26 X 10!° dynes/cm’, | 
14.37 K 10'° dynes/cm’, | 


= / 9 - ri / » f 
7.09 X 10" dynes/cm?, — ¢i4 = 16.87 XK 10'* dynes/cm’, | 


€33 = 105.67 X 10'° dynes/cm?, tie 


wn 


Ca = 
p = 2.654 grams/cm’, } 


while the relations between applied electric fields and strains are: 


j 


= d,E,, e::=- - di4Ey, | 


; i 
Cyy = —dyE,, ery = — 2dyEy, | 


Cre 


are me 
dy = 6.45 X 10° 8. units, 


(23) 
1.45 X 10-°C.G.S. units. | 


— dis 


(1) In X-cut quartz plate designated by X in Fig. 1, if the boundary 
surfaces be expressed by x =0 and x =a, (8) reduces to 


C1— Cc 0 f) 


ul 
So 


0 $ (C4) — Ci) — C Ci4 (24) 
0 Ci4 Coa = 6 


2 W. Voigt, Lehrbuch der Kristallphysik, 754 and 789 (1928). 
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or 
C, = ¢n = 85.45 X 10! dynes/cm?, 
C2, C3 = 4(C11 — Ci2) + ay + [}4(cu — C99) = deus}? + C14]! . (25) 
co = 67.22 X 10! dynes/cm?, 
c3 = 28.98 XK 10'° dynes/cm?, 
and the ratio 1:m:n corresponding to these values: 
Liimyin, = 1:0:0, 
lotmeinNe = Otce — Cuicuy = O: — 0.6021, > (26) 
l3:m3:N3 = O23 — Cagicig = 021.6721. ) 
The normal modes of vibration are given by the following equations: 


— = A cos (qrx/a)-e'”, (27) 


where £;=u; &= —0.60 v+w; and &;=1.67 v+w; while the electric field in 
the direction x cannot excite the displacement & and &; as is seen from (23), 
so that the period of vibration is 


Pp ‘24 = (q, 2a)(¢1;/p)''? = qg/a X 0.2837 & 10° cycles, Sec... (28) 
and the mode of vibration is 
u = A cos (qrx/a)e'”' (29) 


that is, the pure longitudinal vibration along the normal to the surfaces. 

As an experimental example, a quartz oscillating crystal 0.0922 «2.610 
X< 2.684 cm* (the normal to the major surfaces inclines 7’ the axis x) had the 
fundamental frequency of thickness vibration 3.1010® cycles/sec., from 
from which we get 


p/2m = (1/a) + 0.0922 XK 3.10 K 10° = (1/a) XK 0.286 XK 10° cycles/sec., (30) 


which is seen to be very near upon the theoretical result (28). 
(2) In Y-cut quartz plate designated by Y in Fig. 1, if the boundary 
surfaces be expressed by y=0 and y=a, (17) reduces to 


| (Ci ta) = € 0 0 | 
| 0 Cu-¢ =O ey | =O (31) 
0 == €ne Cegu C | 


or 

(= car — C12) = 39.10 X 10'° dynes, cm?, 
Gay a C44) + [2 (cx = C44)" + C142]! ’ 
ce = 93.31 X 10° dynes/cm?, 


cz; = 49.23 & 10"° dynes/cm?, 


C2; C3 = 


and the ratio /:m:n corresponding to these values: 


— 
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liimyin, = 1:0:0 
lotmeinme = O24, — Cotcig = 0:2.15:1 (33) 
lg: m3inzg = O24, — C3iCig = O: — 0.4721 
The normal modes of vibration are given by 
§ = A cos (qry/a)-e'” (34) 
where £,=u, &2=2.157+w, and &; = —0.47 v+w, of which £ and £3 cannot be 


excited piezoelectrically by the electric field in the direction of y axis as is 
clear from (23), so that the displacement is given by 


u = A cos (qry/a)-e'" (35) 


which shows that the vibration is of a pure shear as is pointed out by Pro- 
fessor Cady.* 
The frequency of vibration is 


p q (eu — C12\'” q 

a SE saree yee = x 0.1919 K 10° cycles/sec., (36) 
2x 2a 2p a 

while in the experiment with a sample 0.1909 x 3.200 X 3.304 cm* (the normal 

to the major surfaces inclines 6’ to the axis y) the fundamental frequency of 

thickness vibration is 1.02 10° cycles/sec., that is 


p/2x = 1/a + 0.1909 K 1.02 K 10° = 1/a K 0.195 XK 10* cycles/sec., (37) 


which is also seen to be very near upon the theoretical result (36). 

(3) In quartz plate cut parallel to a surface of positive rhombohedron 
designated by 7 in Fig. 1 (say R-cut plate‘), the vibration is observed to be 
extremely vigorous and the third harmonic vibration (¢=3) can often be 
sustained by Pierce circuit,® in which an oscillating crystal plate is placed be- 
tween grid and filament of a three-electrode thermionic tube. 

To obtain periods and normal modes of thickness vibration, we may first 
transfer the coordinate axes. Let x’, y’, 2’, be so chosen that the direction 
cosines of them referred to the original axes be expressed by the orthogonal 
scheme 





x y Z 
x’ 11 0 0 
: (38) 
y’| 0 cosé —sin#é 
2’ | 0 siné@ cosé 
where 
6 = arctan (3!/2/2.200), (39) 


3 W. G. Cady, A shear mode of crystal vibration, Phys. Rev. 29, 617A (1927). 

* 1. Koga, R-cut quartz oscillating crystal plate and the harmonic vibration, Supple- 
mentary Issue, Jour. I.E.E. (Japan) to be issued in April (1932). 

°G. W. Pierce, Piezoelectric crystal resonators and crystal oscillators applied to the pre- 
cision calibration of wave meters, Proc. A.A.A.S. 59, 81-106 (1923). 
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then the axis y’ becomes normal to the surface r as shown in Fig. 2. If the 
elastic constants referred to the new axes be denoted by c’r., values of c’s in 
(14) are the roots of 


Coe’ =< 0 0 


0 C22, — € co’ |=0 (40) 


or 
C= Cee = C44 sin? 6 ~ R(ciy — C392) cos” 6 =— Css sin 26 
= 62.39 K 10'° dynes/cm’, f (41) 


' 
a , - , Or - / ‘i /\2 _ 49 2 
a(Coe +a) [3 (coe — Cay )® + Co, |! J 


-_ 


Ce, C3 = 
and the ratio /: m:n corresponding to these values: 
liim in, = 1:0:0 ) 
lo:mo:ne = 0:64’ — Corces’ tf (42) 
lsimsins = Otcas’ — C3iCos’ 
The normal modes of vibration are given by 
— = A cos (qry’ a)-e'? (43) 


where £, =u’, £2 =(C44’—C2) v' +004'w’, £3 = (Cag — C3) v' +C24’w’. On the other 
hand, the relations (24) transform into the following relations: 


C2’: = di,E,, ) 


Cy'y’ = — cos 6(d;, cos8 + diy sin #)E,, 
€22 = — sin 6(d;; sin @ + d,4 cos O@)E,-, 

: ' (44) 
€yz’ = (dy, sin 26 + d,,)E,y, 


€xz' = — (discos 6 — 2d; sin 6)(E, cos 6+ E, sin 8), 
Czy = (dissin @ — 2d,, cos 6)(E,y cos 6 + E, sin 6), } 


from which we see that we cannot excite the displacements £ and &; by the 
electric field along the axis y’, consequently the displacement is given by 


u’ = A cos (gry’/a)-e'" (45) 


which shows that the mode of vibration is of a pure shear. The frequency of 
vibration is 
p le = (g/2a)(cee'/p)''? = g/a X 0.2424 X 10° cycles/sec., (46) 


while in the experiment with a sample 0.1002 2.592 X 2.762 cm? (the normal 
to the major surfaces inclines within 1’ to the axis y’) the fundamental fre- 
quency of thickness vibration is 2.47 X 10° cycles/sec., that is 


p/ 2x = 1,4 X 0.1002 X 2.47 K 10° = 1/a X 0.247 X 10® cycles/sec., (47) 


which is also seen to be very near the theoretical result (46). 
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There have been published several papers® on the excitation of harmonic 
vibrations in a resonator, but as far as we are aware, very little work has 
been done on the harmonics of the thickness vibration of a crystal. The 
present crystal plate can often be sustained in its first and third harmonic 
vibrations separately by the Pierce circuit above mentioned. Fig. 3 shows the 
characteristic features of the oscillator. The higher frequency of the vibra- 
tions is easily observed by a heterodyne frequency-meter to be just three 
times that of the lower. 








Frequency , 110 
0.15910" cycles - 


























is 
4 Variable Air Condenser 
| =) = (?) 1500 ppF in 2500 div. t 
| at “Weston Thermoammeter ' {T 
d Hire Model 507 500 mA | 
Quartz Oscillating t ~ 
Crystal Plote D.C.Milliammeter " 
| | O45S= 2.67228) em? | 
—C in div. 
1100 hgh, 
Fig. 3. 


Examples of dimensions which are suited for the third harmonic vibra- 
tions are: 





Dimensions (cm*) Approx. frequency of vibrations (cycles/sec. ) 
0.093 K 2.56 XK 2.57 2.66 X 10° and 2.66 X 10° X 3 
0.114 K 2.09 K 2.70 2.17 X 10° and 2.17 X 10° x 3 , (48) 
0.155 K 2.47 XK 2.81 1.59 X 10® and 1.59 X 10° K 3 
0.185 K 2.77 K 2.98 1.34 & 10° and 1.34 X 10® X 3. ) 


(Inclinations of the major surfaces to the x’z’ plane are generally not greater 
than 10’). 

(4) If the boundary surfaces are parallel to surfaces of a negative rhom- 
bohedron designated by r’ in Fig. 1, normal modes and frequencies of vibra- 
tion can be studied as in the previous case. The only difference is the reversal 


of the sign of @ in Fig. 2. Thus 
Cos. = C44 Sin? 6 + 3(Cu — C12) cos? 8 + ci4 sin 20 (49) 
= 29.59 & 10! dynes/cm?, 


and 
p/ 2m = (q/2a)(ces’/p)'!? = g/a X 0.1670 X 10° cycles/sec. (50) 


8’ E. Giebe and A. Scheibe, Sichtbarmachung von hochfrequenten Longitudinalschwin- 
gungen piezo-elektrischer Kristallstaebe, Zeits. f. Physik 33, 335 (1925). 
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The displacement 
u’ = A cos (gry’/a) -e'”' (51) 


gives a pure shear vibration. 

From the experiment with a sample 0.0765 K 2.76 X3.19 cm* (the normal 
to the major surfaces inclines 5’ to the axis y’) the frequency of vibration is 
2.20 K 10° cycles /sec., that is 


p/ 2x = 1/a X 0.0765 X 2.20 X 10° = 1/a XK 0.168 XK 10° cycles/sec., (52) 


which is very near the theoretical result (50). 

It is seen above that the frequency of vibration of a piezoelectric oscillat- 
ing crystal is always a little higher than that calculated from the elastic con- 
stants. This small difference must be, to a certain extent, due to the finiteness 
of the boundary surfaces of plate, but it is not overlooked that the frequency 
of pure mechanical vibration is different from that as a piezoelectric oscillat- 
ing crystal. In fact, when the crystal plate is excited to vibrate with its 
natural vibration, the equivalent elastic constants become slightly greater 
than that which is determined elastostatically. For example, by a strain e,,, 
besides the stress cye,,, the polarization of €ye,, is produced, €,, being the 
piezoelectric constant, so that the crystal, being acted upon by the electric 
field 47r€,,e,,/K, where K is the dielectric constant, has a stress 47€);e,,/K, 
which adds to the above-mentioned stress cy,e,,. Consequently the resultant 
adiabatic elastic constant c¢,,° becomes 

dire,’ 4dr X (4.77 X 10*)? 
Cn = On +o — = ¢) 41 + — P= (1 + : 0.0074). (53) 
K ("85.45 x 10% x 4.5f 

In other words, to replace a quartz oscillating crystal by an equivalent 
electrical network,’ we must, as shown in Fig. 4, consider C;, which represents 
the capacity of the quartz as a mere dielectric, parallel to the series resonance 


-» 








Fig. 4. 


circuit Lo, Co, Ro representing the mechanical vibrating system. If such is the 
case, we can easily understand that the electrically measured natural fre- 
quency of vibration determined by Lo, Co, Ci, is always higher than that de- 
termined by Lo, Cy. The ratio Co: C,, is in general, not greater than 1/100. 


7 D. W. Dye, Piezoelectric quartz resonator and equivalent electrical circuit, Proc. Phys. 
Soc. London 38, 399-457 (1926). 
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IV. THICKNESS VIBRATION OF TOURMALINE PLATE 


Tourmaline is not yet widely used as a piezoelectric oscillating crystal, 
and we have no occasion to investigate by fine crystal, but we can foresee 
the normal modes and frequencies of vibration in the first approximation. 

The adiabatic elastic constants referred to the coordinate axes shown in 
Fig. 5 are composed of the same scheme with (22), the numerical values being 


ci = 270.17 X 10" dynes/cm?, C12 69.06 & 10!° dynes/cm?, 


C33 = 160.69 X 10! dynes/cm?, (3 = 8.83 X 10'* dynes/cm’, , 
Cy = 66.71 XK 10" dynes/cem?, — ¢44 = 7.75 X 10! dynes/cm?, = 
p = 3.100 grams/cm* } 
while the relations between electric fields and strains are: 
Crz = — dak, + dziE2, ey: = di;Ey, 
Cyy = dwk, + d3E., €zze = 4;E:, 
€.2 = d33F., Cry = — 2dekz, » (55) 


— dx =0.69 K 10-5 C.G.S. units, d3; = 0.74 K 10°°C.G.S. units 
d33 = 5.78 X 10°*C.G.S. units, d);; = 11.04 * 10-° C.G.S. units | 


(1) If the boundary surfaces of plate be z=0 and z=a, (20) reduces to 


C4ag — C 0 0 
0 ta~e 0 = 0) (56) 
0 0 C33 — C 
or 
C, = Co = C44 = 66.71 X 10! dynes/cm’, 1 - 
C3 = C33 = 160.69 X 10'° dynes/cm?, f ef) 
and the ratio /:m:n corresponding to these values: 
1,2 m,:n, = 1:0:0 ) 
lo? met Me = 0:1:0 + (58) 
ls 2 m3: M3 = 0:0:1 
The normal modes of vibration are given by the following equations: 
£; = u = A cos (qrx/a)-e'", 
fo = v = A cos (qrx/a)-e'", , (59) 
&; = w = A cos (qrx/a)-e'", 


of which w can only be excited piezoelectrically by the electric field along the 
axis Z, as is seen from (55). Accordingly the vibration is of a pure longitudinal. 
The corresponding frequency of vibration is 


p/2a = (g/2a)(c33/p)'/* = (q/a) X 0.3600 XK 10° cycles/sec. (60) 
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(2) If the boundary surfaces of plate be x=0 and x =a, (8) reduces to 
C1— Cc 0 0 | 
0 $(¢11. — Ci2) — € C14 = () (61) 
UV Ci Cam Cc | 
or 


CQ, = Cy = 270.17 x 10!° dynes, cm?, 


an mae  - 1 - te — 2 
Ce, C3 = (Cu — C12) + 3€ag + [$3 (cx — (33) - 20444 + C14 }} 


(62) 
ce = 102.25 X 10'° dynes/cm?, 
cz; = 65.03 XK 10'° dynes/cm’, 
and the ratio /:m:n corresponding to these values: 
litmyin, = 1:0:0 } 
lo: meine = O2Ce — Cygidig = O: — 4.5921 + (63) 
l3:m3:m3 = O2¢3 — Cagicig = 0:0.22:1 ) 
The normal modes of vibration are given by the following equations: 
— = A cos (qrx/a)-e'” (64) 
where £;=u, &2= —4.59 v+w, &£;=0.22 v+w, while the electric field in the 


direction x cannot excite the displacement £, and can hardly excite £2 as is 
seen from (56), and the periods of £ and &; are: 


q/a X 0.2872 X 10° cycles/sec., (65) 
p/2a = (g/2a)(c3/p)'/* = g/a X 0.2290 XK 10° cycles/sec. (66) 


p/2a = (q/2a)(c2/p)'!? 
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Resistance-Temperature Relation of Beryllium 
Oxide and Zinc Oxide 


By Joun A. OSTEEN 
Department of Physics, Lehigh University 


(Received June 13, 1932) 


Bidwell showed that the resistance-temperature law for certain variable con- 
ductors (oxides, etc.) was of the form p = Ae®/#Tte7, Plotted in the form (1/p)(dp/dT) 
= (/RT?—a he obtained for Fe,Q; straight lines with a break near the recalescence 
point, the two lines having the same y intercept but different slopes. This law was 
also followed by metallic germanium. The present paper extends this work to zinc 
oxide and beryllium oxide. Beryllium oxide yields a straight line showing a transition 
at 750°C. Zinc oxide yields two straight lines with a transition in the interval 250°C— 
500°C. The two zinc oxide lines have different slopes but the same y intercepts. With 
successive heatings to 1050°C there occurs continued decreases in the slope of the line 
for the 500°C-1000°C range but no change in the y intercept. The line for the range 
0°C-250°C changes on successive heatings to 1050°C and finally stabilizes to a value 
of slope and intercept which repeats on succeeding runs. On the suggesion that the 
indicated transformation was due to impurity, new material of special purity pre- 
pared by the New Jersey Zinc Company was studied. The behavior of this material 
was found to agree closely with that of the earlier specimen but gave a more clean 
cut transformation. 


INTRODUCTION 


ONIGSBERGER' on the basis of the electron theory of Drude derived, 
by making use of certain empirical relations, the expression 


p = A(1i + at + Bi*)e@/*7 (1) 


for the resistance of variable conductors. (p) is the specific resistance (used 
in this work as total resistance); (¢) the temperature on the centigrade scale; 
(7) the temperature on the absolute scale; (R) the gas constant (approxi- 
mately 2); (Q) the heat of dissociation of the electrons, i.e., the heat required 
to liberate a gram-atom of electrons from the atoms to which they are sup- 
posed to be bound; (e) the Naperian base; (A), (@), and (8), constants. 

In the case of Fe,O; Bidwell? showed that the resistance-temperature 
relation could be written in the form 


- 


p = Ae@/RT+aT , (2) 


Eq. (2) differs from (1) in that the exponential e*7 is substituted for (1+ at 
+6t?). This proved to be empirically satisfactory and more convenient 
than (1) for graphical testing. 

By expressing Eq. (2) in the form 


logiop = 0.435(0/RT + aT) + logio A (3) 
and differentiating, we have 


1 Kénigsberger, Ann. d. Physik 32, 179 (1910). 
2 Bidwell, Phys. Rev. (2) 8, 12 (1916). 
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[— (1/p)(dp/dT)| = Q/RT? — a. (4) 


If a straight line is obtained by plotting the left hand member of Eq. (4) as 
ordinate against 1/7* as abscissa, the equation is verified. 

Bidwell obtained for Fe.O; straight lines with a break in the 600°C to 
800°C region. This break he relates to the recalescence point in iron which 
occurs in this region. The present paper offers some results for the resist- 
ance-temperature relations for specimens of BeO and ZnO as plotted by the 
above method. In addition, for comparison, logarithmic plots, similar to 
those usually made for variable conductors are also included.* 


PREPARATION OF SPECIMENS AND METHODS EMPLOYED 
The materials used were C.P. BeO prepared at Cornell University; Bakers 
C.P. ZnO and S.P. ZnO prepared by the New Jersey Zinc Company. 
Data were obtained on compressed and baked specimens prepared as ex- 
plained by Bidwell for iron oxide.” 
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DISCUSSION OF DATA FOR BEO 


Fig. 1 shows for two runs on the specimen of BeO resistance plotted di- 
rectly against temperature. The resistance value at any point is obtained by 
multiplying the reading taken on the resistance-axis by the indicated power 
of 10 of the curve from which the point is taken. The heating interval between 
runs (1) and (2) did not affect the absolute values as is seen in the coincidence 
of the curves. From this curve, tangents were drawn as shown and with the 
data thus obtained Eq. (4) was plotted. This gave a straight line from 275°C 
to 750°C. The value of Q (17,800) from this line is of the order of values re- 


* Spanner, Ann. d. Physik 75, 609 (1924), 
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ported by Kénigsberger for CuO, SiO., ThO and CeO, as well as near the 
value 17,400 obtained by Bidwell for Fe.O3. 

The noticeable lack of smoothness in the resistance-temperature curve 
between 750°C and 800°C is followed, in the plot of Eq. (4) by a transition 
the last three points of which indicate another possible straight line with re- 
sistance intercept of same value as the 275°C to 750°C line. 
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Since the equation is satisfied, the conductivity appears, an the basis of 
the derivation of the equation tested, to be electronic in the 275°C to 750°C 
range. However, some polarization which occurred in the runs indicates that 
electrolysis is also involved in the conductivity. 
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Fig. 2 shows the logarithmic plot for specimen one. Here as in the plot of 
Eq. (4) a straight line is obtained from 257°C to 750°C with points erratic 
beyond 750°C. 

The plot of —(1/p)(dp/d7) against 1/7? for a second specimen is shown on 
Fig. 1. This specimen, baked in a platinum furnace, gave a constant value of 
resistance for 1500°C after six hours heating. The Eq. (4) in this case is fol- 
lowed for the 1160°C to 1360°C range. 

A second run on specimen 2 for a 980°C to 1530°C range gave points that 
were rather erratic for the test of Eq. (4). Fig. 3 shows the results for this 
range for log conductivity plotted against 1/7. 

As a check on the break occurring at 1227°C two runs of thermoelectric 
power plotted against average temperature is shown in Fig. 4. The 1227°C 
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Fig. 5. 
break correlates very well with the 1200°C thermoelectric power break. 
Fig. 3 also shows the 1160°C to 1360°C run plotted logarithmically. 


DiIsCUSSION OF DATA FOR ZNO 


Fig. 5 is the third of the resistance-temperature runs for a specimen of 
C.P. ZnO. The circles indicate the points obtained with decreasing tempera- 
ture, the crosses with increasing temperature. This partial run from room 
temperature to about 400°C was made in order to check the peculiar curva- 
ture in the 250°C to 300°C range. Each of the three runs showed definitely 
this peculiarity. 

A preliminary run made after the specimen had been baked for thirty 
minutes in the iron container in which it had been compressed gave quite 
erratic results in particular throughout the region 550°C to 250°C. The speci- 
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men was therefore baked for about four hours at approximately 1000°C. 
After this, runs one, two, and three were made which gave such curves as 
shown in Fig. 5. The heating interval necessary between runs evidently 
changed the condition of the specimen since these three curves differed in 
absolute values of resistance. 
From the curves of the three runs tangents were drawn and the data 
plotted according to Eq. (4), (Fig. 6). 
The three plots give straight lines from 550°C to 950°C and from 25°C 
to 175°C with a transition from 175°C to 550°C. The two lines have different 
slopes, but the same resistance intercepts. The plots of these runs show that 
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with successive heatings to 1050°C there occurs continued decreases in the 
slope of the line for the 500°C to 1000°C range but no change in the resistance 
intercept. The line for the 25°C to 175°C range repeated for values of slope 
and intercept after run one. 

On the suggestion of H. M. Cyr of the New Jersey Zinc Company that 
the transition indicated in the 175°C to 550°C range might be due to im- 
purity, a new sample of S. P. ZnO supplied by the New Jersey Zinc Company 
was studied. 

Fig. 7 shows both the temperature-resistance curve and the plot for 


[— (1/p)(dp/dT)] = O/RT? — a. 


The break in the 250°C to 300°C range does not clearly appear in the tem- 
perature-resistance curve as it did in the C.P. specimen. However, the plot 
of [—(1/p)(dp/dT)] against 1/7? indicates an irregularity in the 200°C to 
550°C range similar to that obtained in the C.P. specimen. 
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This irregularity in the 200°C to 550°C range seems to be related to the 
growth of particles since below 500°C some change apparently occurs as 
rubber tests, which are very sensitive to particle-size, show measurable dif- 
ferences if the oxides have been heated to 300°C. There is also an apparent 
electronic change in this region which is indicated by the yellow color of 
zinc oxide. Experiments by the New Jersey Zinc Company have shown that 
at 350°C a yellow tint is just noticeable to the eye. On further heating, the 
color gradually brightens and becomes a deep pure yellow at 500°C. After 
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500°C a more stable condition is reached as indicated by the fact that Eq. 
(4) follows a straight line above 550°C. 

In connection with the yellow tint mentioned above, it was found that 
the S.P. ZnO subjected to the 10,000 pounds per square inch pressure be- 
came lightly yellow even for room temperature. On heating and then cooling 
to room temperature the material became white again. 

The above runs were made in a nichrome wound furnace. Fig. 8 shows the 
temperature-resistance run for the S.P. specimen of ZnO when run in a 
platinum-wound furnace. The temperature about the specimen was raised 
to approximately 1500°C, where the resistance measured about one ohm. On 
decreasing the temperature the resistance increased very slowly reaching a 
value of 1.4 ohms at 470°C. Thereafter the resistance increased rapidly 














RESISTANCE-TEMPERATURE RELATION 89 


enough to give the curve of Fig. 8. From this curve tangents were drawn and 
the plot of |—(1/p)(dp/dT)]=Q/RT?—a gave the lower straight line shown 
in Fig. 7, which has a different slope but same intercept as the 925°C to 
550°C and 200°C to 25°C lines. 

The specimen during the above run shrank and formed quite noticeable 
crystals on the surface. The color was green. The specimen also showed the 
property of rectification. This property is the same as that found by the New 
Jersey Zinc Company for clinkers of zinc oxide taken from their furnaces 
and used in earlier radio work for rectifiers. 

Fig. 9 is the logarithmic plot for both the C.P. and S.P. specimens. The 
transitions in these correlate very well with those of the Eq. (4) plot. 

The author is indebted to Professor C. C. Bidwell for supervision in this 
work, and to Mr. H. M. Cyr of the New Jersey Zinc Company for the S.P. 
zinc oxide. 
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Recent published results of studies on the flow of liquids through porous media 
under the action of gravity have shown wide disagreement in both the formulation 
and interpretation of the problem. A new attack on the problem has been carried out 
and has led to unambiguous answers to the questions of interest. Experiments on a 
radial sector of sand show that the old Dupuit formula of 1863 stating that the fluid 
outflow is proportional to the square of the differences in the fluid heights in the sand 
is exact within experimental error provided the fluid heights are replaced by fluid 
heads as measured at the sand bottom. Both the pressure distribution formula and the 
integrated expression for the fluid outflow are verified in detail. The cases where an 
added pressure head is superposed upon the gravity flow, heretofore not mentioned in 
the literature, behaves both with regard to its pressure distribution and fluid outflow 
as a direct superposition of the simple gravity flow and pure radial flow systems. A 
semi-quantitative treatment is given for the flow in the capillary zone above the main 
fluid body. It is necessary under certain experimental conditions to correct for the 
added flow carried by the capillary layer in addition to that induced by the gravity 
drive. In most practical cases the capillary flow may be ignored. A theoretical dis- 
cussion is given of the conditions at the free surface of a gravity flow system and it is 
shown that nowhere on it can the slope exceed 45°. The results on the radial flow ex- 
periments are also generalized to give a method for treating general drives. 


I. PURPOSE OF INVESTIGATION AND CRITIQUE OF PREVIOUS INVESTIGATIONS 


T HAS very frequently been assumed that the effect of gravity is an en- 

tirely negligible factor in the production of oil from sands. It is true that 
during the early period of production while the sand pressures are fairly high 
and the pressure at the well is maintained relatively high, the contribution of 
gravity flow to the total is negligible or entirely absent. However, when the 
reservoir is depleted of its gas and the oil has lost its expulsive powers due to 
depletion of dissolved gases, one cannot arbitrarily ignore the effect of gravity 
on production rates. To derive a definite answer to this problem was the aim 
of the following investigation since, as will be seen presently, the results 
presented in the available literature on the subject can hardly be considered 
as satisfactory. 

Although the problem of gravity flow does not seem to have been con- 
sidered from the point of view of its relation to oil production it has been 
given considerable attention hy investigators interested in artesian flow and 
irrigation problems. Nevertheless, a review of the literature shows that 
neither froman experimental nor froma theoretical point of view has a definite 
solution to the problem been obtained. 

Thus, as long ago as 1863 Dupuit! made the assumption that for small 
inclinations of the free surface the stream lines can be taken as horizontal and 


1]. Dupuit. Etudes Theoriques et Pratiques sur le Mouvement des Eaux. 
i I ] | 
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are associated with velocities independent of the depth, being equal to that 
at the free surface (shell flow cf. §6). In particular he derived a formula 
applying to radial gravity flow. This was then generalized by Forchheimer’ 
so as to give a general equation for the free surface of any gravity flow system 
(cf. §7 Eq. (20)). Although the assumptions underlying these equations have 
but little validity in regions surrounding the surfaces of outflow of the sys- 
tems (cf. §6), they have been used indiscriminately by Forchheimer and 
others in solving special problems, apparently without any attempt to test 
empirically the predictions they found. 

Perhaps the most satisfactory study of any problem of gravity flow in 
porous media is the elegant treatment by Hopf and Trefftz* of the problem 
of the partial drainage by a ditch of the regional gravity flow in a water sand 
with an inclined impermeable bed. This problem, idealized as a two-dimen- 
sional one, was solved exactly by these authors by using the method of conju- 
gate functions as applied to the rigorously valid Laplace’s equation in the 
interior of the sand (cf. §6). Unfortunately, however, this method cannot be 
extended to three dimensional cases such as the radial flow problem studied 
in the present investigation. 

Another rigorous attack on the problem of gravity flow, that is, one using 
solutions of the potential equation (cf. Eq. (13a) ), was made even earlier by 
J. Kozeny* with special reference to the case of radial flow. Unfortunately, 
however, Kozeny neglected entirely the discontinuity between the fluid 
height in the sand at the well face and the fluid level in the well. Thus he 
required of his uppermost streamline, which represented the total flux 
through the system, that it leave the sand body at the fluid level in the well. 
Further, Kozeny made no attempt to make his uppermost or free surface 
streamline correspond to a real free surface over which the potential function 
at any point is directly proportional to the vertical coordinate of that point. 
Then, too, Kozeny’s result that the flow through the system is proportional to 
the well radius depends on the assumption that the well radius is large as 
compared to the maximum fluid head in the sand. Just the converse, how- 
ever, corresponds to most cases of practical interest and when it is taken into 
account even Kozeny’s formulae give an inverse logarithmic variation with 
the well radius, as is found in the present investigation (cf. §4, Eq. (4)). 

It therefore appears that although Kozeny’s method of attack is quite 
rigorous, the assumptions made in its development must necessarily vitiate 
the validity of his conclusions. Thus it is not surprising that he gets the 
peculiar result that the flow through a radial flow system is a maximum when 
the fluid head within the well is half of that at the external boundary. Like- 
wise, Kozeny’s result that the flow would vanish with vanishing fluid head 
in the well is clearly due to the neglect of the discontinuity at the well. 

The apparent check of Kozeny’s results by Kozeny® himself or the more 


2 Ph. Forchheimer, Zeitschr. d. Arch. u. Ing. Ver. Hannover (1886). 
* Hopf and Trefftz; Zeitschr. f. angew. Math. u. Mech. 1, 290 (1928). 
4 J. Kozeny. Wasserkraft u. Wasserwirtschaft, 22, 120 (1927). 

5 J. Kozeny. Wasserkraft u. Wasserwirtschaft 22, 146 (1927). 
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extensive experiments by R. Ehrenberger,® however, do not seem conclusive. 
Thus the point in Kozeny’s empirical data which he takes as the indication 
of the maximal Q for the fluid well height, being half of that at the back of 
the sand, might well be taken as being due to experimental error since it can- 
not lie upon any smooth curve drawn to represent a major part of the data. 
Since for comparison the Dupuit equation (cf. Eq. (18)) was forced to pass 
through this point it cannot be expected that it would check with the other 
points that do lie on a smooth parabolic curve. 

Ehrenberger’s experiments, though beautifully performed and exhaustive 
in many respects, do not seem to give a satisfactory answer to the problem. 
In his desire to correlate his results with Kozeny’s theory, Ehrenberger was 
unfortunately misled in his method of analysis. Thus when he finds that 
Kozeny’s formula does not work he sets up an empirical correction factor 
which evidently has no physical significance. If he had replotted his Fig. 20 as 
QO vs. AZ*?, as suggested by the Dupuit formula and as done in the following 
discussion, rather than as Q vs. Z,, as suggested by Kozeny’s formula, Ehren- 
berger would have found a straight line confirming the present investigation 
and would not have had to resort to the empirical correction factor. 

It will be clear, therefore, that in spite of the many studies already 
carried out on the problem of gravity flow in liquids a satisfactory description 
of the problem had not been attained. The investigation to be described be- 
low was undertaken in an attempt to clear up the situation and to derive a 
conclusive answer with regard to as many features of the problem as possible. 


I]. EXPERIMENTAL EQUIPMENT AND PROCEDURE 


The principle of the experiments was essentially that of constructing a 
section of a model well flowing under the action of gravity and observing the 
production rate under various conditions of flow. It should be mentioned at 
this time that while the original equipment was designed for gravity flow only, 
it was found imperative that permeability’ measurements be made on the 
sand in place in the model, and a second model was therefore constructed. 
The later model was pressure tight so that from strictly radial flow measure- 
ments the actual permeability of the sand could be accurately measured. In 
addition this model permitted an experimental check on the composite flow 
formula deduced from a knowledge of the radial and gravity systems. 

The initial gravity flow setup was that indicated in Fig. 1. While the same 
results were obtained from the second and more complete apparatus, a dis- 
cussion of the former is included because of the more complete data obtained 
on pressure distributions together with some observations on the behavior of 
the fluid in the capillary zone. 

In Fig. 1, the sector of 15° represents 1/24 of a complete radial flow 
system. The 152.4 cm length of the sector corresponds to the 152.4 cm closest 
to a well of radius 5.95 cm which is the distance at which the screen simulat- 


® R. Ehrenberger; Zeitschr. des Oster. Ing. u. Arch. Ver. Heft 9/10, 11/12, 13/14 (1928). 

7 By the term permeability is meant the fluid outflow through a unit cross section under 
the action of a unit pressure gradient. It is a constant of the sand and liquid system and in- 
versely proportional to the viscosity of the liquid. 
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ing the well face was placed. The model was made covering only the im- 
mediate neighborhood of the well in order to be able to study in particular 
the detailed effects of the convergence near the well and the fluid heights in 
the well, on the free liquid surface within the sand and on production rates. 
It is clear that at great distances from the well the free liquid surface will be 
nearly horizontal and the flow will be essentially of a radial horizontal char- 
acter. This portion, representing the region remote from a well is simulated 
by the sand-free tank used as a reservoir in the model. Thus except for the 
height of the sand (40.5 cm), necessarily small in the laboratory, the setup 
represents closely a full scale model. 

The sand reservoir was separated from the sand-free chamber in the back 
by a fine screen supported against a brass plate perforated uniformly by 
1/2’’ holes with staggered 5/8’’ separations. This insured a uniform distribu- 
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Fig. 1. Diagram of apparatus for radial gravity flow experiments. 
tion of liquid against the far end of the sector. The sand-free tank continued 
12 cm beyond the screen and plate. The glass sides were sealed to the brass 
frame with beeswax-rosin mixture. The sand used was a fairly uniform 
Ottawa sand of about 30/inch mesh. From the preliminary tests with an 
earlier setup it was learned that unless special precautions were taken in fill- 
ing the sector with the sand and water, air bubbles would be trapped in the 
sand pores and destroy the homogeneity of the pore openings for the flow of 
the liquid. The procedure finally adopted in filling the sector was that of 
thoroughly mixing and stirring the sand by hand while keeping the water 
level in the tank just above the sand. Each layer of sand was freed of air be- 
fore more was added. From the uniformity of appearance of the sand as seen 
through the glass sides.one could tell if the sand were really free of trapped air. 
From holes in the bottom of the sector placed along the central radius, 
twenty small manometer tubes were led and mounted next to one of the glass 
sides at points corresponding to their entries in the bottom. The ten tubes 
nearest the well were placed at one half the separations of those nearer the far 
end since the pressure gradients change most rapidly near the well and it is 
here that the greatest detail was desired. These tubes, of course, indicate 
hydrostatic pressures at the bottom of the sand. It was hoped, and in fact 
turned out to be true, that these effective pressures are all that are required to 
predict completely the production from the well. 
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The outlet from the bottom of the well was furnished with a stop-cock 
so as to control the height of fluid in the well and hence the production 
rate. A pump was connected to the outlet collecting vessel to pump the 
outflow water back into the inflow tank, thus maintaining a steady condition 
of flow for given levels in the tank and well. This circulation eliminated the 
accumulation of air bubbles in the sand which would occur if fresh water were 
continually supplied to the system. 

The experiments consisted essentially in noting the rate of outflow of 
water from the sand for given heights of the water level in the tank at the 
back of the sector and in the well at the apex. Readings were also taken of 
the heights in the various intermediate manometer tubes. The experiments 
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Fig. 2. Diagram of apparatus for composite (radial drive plus 
gravity) flow experiments. 


were so conducted that the water level was never appreciably raised into an 
upper zone containing air. This precaution was to confine the flow to air-free 
sand and when it was necessary to increase fluid heights the tank was 
emptied of sand and properly refilled. 

While good use was made of the pressure distributions obtained from this 
setup it was found impossible to obtain effective permeability values for use in 
computing output volumes. It was therefore necessary to construct a new 
pressure tight flow tank in which permeability tests could be made on the 
sand in place. Before going into further detail as to experimental results we 
will describe briefly the second setup used in the final production rate tests 
and in the work on composite flow. 

The general arrangement of the second flow tank is shown in Fig. 2. This 
was again a 15° sector with a well screen at the apex and a sand-free tank sec- 
tion at the back end. The entire sector was constructed of 6.3 mm brass plate 
reinforced with angles to prevent flexure of the sides. All joints were soldered 
except the top plate which was screwed to the top flanges with a rubber 
gasket interposed. Each end of the flow tank was equipped with a stand pipe 
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some 10 cm in diameter and extending 65 cm above the top of the sand sector 
thus enabling the system to be operated as a purely radial flow system by 
maintaining the fluid level at the apex or well end, above the sand level. Ob- 
viously it was necessary to close in the top of the sand body, and to make 
certain that no free space existed between the top of the sand body and the 
top plate of the tank, three glass standpipes 2.5 cm diameter and 65 cm long 
were cemented into the top plate to serve as sand reservoirs. After care- 
fully packing the sector with sand as previously described until completely 
filled, the upper plate was bolted in place with a rubber gasket to seal the 
flange. Water was then circulated continuously with levels at both ends held 
above the sand height. Continuous jarring of the tank served to settle the 
sand and additional sand to fill any voids created, was fed into the sector from 
the sand reservoirs. This process was continued until the permeability 
reached a minimum value under radial flow tests. This minimum perme- 
ability was also checked against the value obtained in a simple glass tube 
linear-flow apparatus where careful packing and elimination of sand-free 
spaces may be observed visually. By this means we were certain that no sand- 
free channels existed in the sector. Preliminary experiments showed the ab- 
solute necessity of this rather elaborate procedure. 

In addition to manometers for measuring the fluid head at the entrance 
and well ends of the sector, five additional manometers were provided to 
measure intermediate pressures. A variable speed gear-pump was used to 
circulate the water continuously at any desired rate and a carefully calibrated 
orifice-plate flow meter was used to measure flow rates. By this means 
absolutely steady state flow conditions could be maintained. 

As before, the experiments consisted in observing flow rates as a function 
of pressure differentials maintained between the boundaries of the flow sys- 
tem with intermediate pressure readings as a further check. The tests in- 
cluded purely radial flows for effective permeability determinations; flows 
in which the input water level was held above that of the sand, while at 
the well end the fluid level was below the sand height, that is, the composite 
flow case; and the usual gravity flow test in which the fluid levels at the both 
ends were below sand height. In order to eliminate the effect of air bubbles 
the same precautions mentioned above were observed, that is, fluid levels in 
the sand were never raised into an upper zone containing air, all changes in 
water level being made downward. 


III. QUALITATIVE DiscUssION OF RESULTS 
Pressure distribution and flow rate data 


The results of the observations from one set of experiments with the 
sector of Fig. 1 are given in Table I. These data give the pressure distribution 
for various liquid heights in the well and the resultant production rates. A 
similar set of observations is given in Table II, the difference being that in 
the latter set the fluid height at the input end of the sector was varied through 
a considerable range. The purpose of showing both sets of data is to show the 
effect of flow in the capillary layer, a feature which will be brought out later on. 
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TABLE I. Flow tank of Fig. 1. 











| Run Run | Run | Run Run | Run 
No.4 | No.5 | No.6 No.7 | No.8 














Gauge | Pres- | 
No. sure | | | 
head | | 
(cm) | | | | 
Well | 29.10 | 27.85 | 23.30 | 20.20 | 15-10 10.90 | 10.60 | 6.10 
1 | 29.45 | 28.15 | 24.30 | 21.80 | 17.70 | 14.70 | 14.40 | 11.80 
2 | 29.60 | 28.40 | 25.25 | 23.30 | 20.20 | 18.00 | 17.70 | 16.00 
3 | 29.75 | 28.55 | 25.85 | 24.30 | 21.80 | 20.10 | 19.75 | 18.50 
4 | 29.90 | 28.65 | 26.35 | 25.00 | 22.90 | 21.50 | 21.15 | 20.10 
5 | 30.10 | 28.85 | 26.85 | 25.70 | 23.90 | 22.80 | 22.35 | 21.45 
6 | 30.10 | 28.85 | 27.10 | 26.20 | 24.70 | 23.75 | 23.25 | 22.45 
7 | 30.25 | 29.00 | 27.50°| 26.70 | 25.40 | 24.65 | 24.15 | 23.40 
8 | 30.35 | 29.20 | 27.95 | 27.20 | 26.15 | 25.45 | 24.95 | 24.30 
9 | 30.50 | 29.20 | 28.10 | 27.50 | 26.60 | 26.00 | 25.50 | 24.85 
10 | 30.45 | 29.20 | 28.20 | 27.70 | 26.90 | 26.40 | 25.85 | 25.25 
11 | 30.45 | 29.20 | 28.50 | 28.15 | 27.55 | 27.20 | 26.60 | 26.10 
12 | 30.60 | 29.30 | 28.80 | 28.60 | 28.10 | 27.85 | 27.20 | 25.75 
13 | 30.60 | 29.40 | 29.00 | 28.90 | 28.55 | 28.40 | 27.70 | 27.35 
14 | 30.70 | 29.45 | 29.20 | 29.20 | 29.00 | 28.90 | 28.20 | 27.85 
15 | 30.80 | 29.55 | 29.50 | 29.50 | 29.40 | 29.30 | 28.65 | 28.40 
16 | 30.90 | 29.70 | 29.75 | 29.85 | 29.85 | 29.90 | 29.20 | 28.90 
17 | 30.90 | 29.70 | 29.85 | 30.00 | 30.10 | 30.20 | 29.45 | 29.20 
18 | 30.90 | 29.65 | 29.95 | 30.15 | 30.35 | 30.45 | 29.75 | 29.55 
19 | 30.95 | 29.70 | 30.10 | 30.30 | 30.60 | 30.65 | 30.00 | 29.85 
20 | 31.05 | 29.75 | 30.25 bw 55 | 30:90 31.00 | 30.40 | 30.25 
Rate 
cc/sec. | 2.75 | 2.68 | 9.05 | 13.15 | 18.12 | 21.60 | 21.60 | 22.30 





In the second set of experiments the top free surface of the liquid in the 
sand was also noted as accurately as possible and recorded. Two sets of these 
readings which are typical of all these sand height data are plotted in Fig. 3. 
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Details of 


As will be seen later, the data already given are sufficient to give a descrip- 





Run | Run 
No. 9 — 
| 

0 0 

9.90 | 10.15 
15.05 | 15.00 
17.80 | 17.50 
19.45 | 19.15 
20.90 | 20.55 
21.95 | 21.50 
22.95 | 22.50 
23.90 | 23.40 
24.50 | 23.95 
5.00 | 24.35 
25.75 | 25.20 
26.50 | 25.85 
27.05 | 26.40 
27.60 | 26.90 
28.10 | 27.45 
28.65 | 27.95 
28.90 | 28.25 
29.35 | 28.60 
29.70 | 28.90 
30.10 | 29.30 
22.70 | 21.3 
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Fig. 3. Fluid levels and fluid heads in sand (gravity flow). 
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TABLE II. Flow tank of Fig. 1. 





Run Run Run | Run Run | Run Run | Run | Run 
No. 11 | No.12 | No.13 | No. 14 No. 15 No. 16 | No.17 | No. 18 No. 19 
Gauge | Pres- | 






































No. sure | 
head | 
(cm) | 
Well 32.50 27.70 22.50 17.00 16.00 15.60 7.00 | 5.80 0 
1 33.50 29.30 25.00 20.60 19.40 18.20 12.50 11.30 5.80 
2 34.20 30.50 26.90 23.40 22.00 20.20 16.20 15.10 7.30 
3 34.80 31.70 28 .40 25.80 24.10 21.80 18.80 17.60 8.30 
4 35.30 32.60 29.80 27.50 25.70 23.10 20.70 19.40 8.90 
5 35.85 33.40 31.00 29.00 27.00 24.10 22.30 20.80 9.50 
6 36.20 34.00 31.90 30.10 28.00 25.00 23.30 22.00 9.90 
7 36.60 34.60 32.70 31.10 28.90 25.70 24.30 | 22.80 10.30 
8 36.80 34.90 33.20 31.70 29.60 | 26.20} 25.10 23.60 10.60 
9 37 .00 35.40 33.90 32.50 30.20 | 26.60 | 25.80 | 24.10 10.90 
10 37.20 35.60 34.20 33.00 30.50 26.90 | 26.10 24.40 10.90 
11 37.60 36.30 35.00 34.10 31.50 27.60 27.10 25.40 11.20 
12 37.80 36.70 35.60 34.80 32.20 28.10 27.80 26.10 11.60 
13 38.10 37.10 36.20 35.50 32.70 28.60 28.30 26.60 11.70 
14 38.30 37.50 36.70 36.10 33.20 29.00 28.90 27.10 11.90 
15 38.40 37.80 37.20 36.80 34.00 29.60 29.60 27.80 12.30 
16 38.70 38.10 37.50 37.20 34.40 29.80 29.90 28.00 12.40 
17 39.10 38.70 38.20 38.00 34.90 30.30 30.50 28.60 12.70 
18 39.10 38.80 38.50 38.30 35.10 30.40 30.80 28.80 12.70 
19 39.30 39.10 38.80 38.70 35.50 30.60 31.00 29.20 12.90 
20 39.30 39 .30 39.10 39.00 35.90 31.00 31.50 29.50 13.10 
Rate 
cce/sec. | 10.70 17.10 22.90 28.10 24.30 18.00 23.90 21.80 5.79 
Radius—apex of system to gauge (n) =r, 
cm cm cm cm 
lw 6.4 rs 30.5 Ti 66.0 Tie 116.8 
rT, 10.2 '6 35.6 T\2 76.2 "17 127.0 
re 15.3 7 40.6 T13 86.4 ris S372 
T3 20.3 rs 45.7 T14 96.5 rig 147 .3 
%%; 25.4 Ty 50.8 Tis 106.7 Toy 156.0 
8 
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tion of the problem as a whole. Some interesting information as to the details 
of the flow in the sand was obtained in the following manner. A narrow tube 
with a capillary opening at the tip was forced down vertically into the sand 
at the glass face and near the back end of the sector. Small quantities of black 
ink were added through this tube so as to maintain.a small steady flow of ink 
from the tip into the sand. This ink formed streamers that progressed along 
the sector towards the well with surprising sharpness. They clearly proved 
that the flow in the sand was of a streamline character and at the same time 
they traced out these streamlines. 

The detailed positions of these streamlines were recorded and then plot- 
ted. The results are shown in Fig. 4. The minor irregularities are undoubtedly 
due to local inhomogeneities in the sand body. The downward slope of the 
lower streamlines at the tank end may be explained by the fact that the 
lowest row of holes in the back plate was about 1 cm or more above the bot- 
tom so that water entering the sand from the supply reservoir must diverge 
downward as indicated by the streamlines in this region. 
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Of more interest are the marked upward slopes of the upper streamlines 
at the reservoir end of the sector. This may at first seem surprising since 
one might at first suppose that the upper part of the wetted sand would be a 
dead capillary layer taking no appreciable part in the flow. When it is recalled 
however that due to the pressure discontinuity across the capillary free sur- 
face there must be a pressure deficiency immediately below it, it becomes clear 
that the capillary layer instead of being dead should actually act like a sy- 
phon and play a significant part in the flow of the liquid. To do so obviously 
requires that upon entering the sand-body from the supply tank, the stream- 
lines in and nezcr the capillary layer slope upwards as was observed. Further 
discussion as to the amount of the flow occurring in the capillary zone will 
be presented later. 
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Fig. 4. Observed streamlines. 


Besides the character of the streamlines near the tank it is of interest to 
note their nature near the well. Although the data are not sufficiently ac- 
curate to yield exact values of the slopes of the streamlines it is fairly clear 
from the curves that only the very upper ones have inclinations with the 
horizontal exceeding 45°. Further, it is to be noted that the capillary free 
surface itself is much flatter near the well than the streamlines below it and 
in fact goes through an inflexion as is indicated more clearly in Fig. 3. These 
qualitative features fit in directly with some analytical results to be given in 
$6. There it is shown that under no conditions in gravity flow in a homogen- 
eous medium (neglecting kinetic energy) can the slope of a free surface ex- 
ceed the value 1, corresponding to an inclination of 45°. In general, the slope 
of an ordinary streamline will have an upper limit lower than 1. Only if the 
pressure distribution is such that the vertical velocity decreases with the 
height may one get inclinations exceeding 45°. Since there is a pressure de- 
ficiency underneath the capillary free surface it is quite possible that as one 
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goes into the capillary layer the effect of the capillary surface will be to de- 
crease the vertical velocities and thus make the slopes exceed 1. Hence the 
inclinations exceeding 45° are even to be expected for such streamlines. 
Further, since the capillary surface itself cannot exceed 45° the streamlines 
immediately below it, near the well, should break away and leave a space 
relatively dead in the immediate neighborhood of the well and just below the 
free capillary surface. This too, corresponds to the observations as is in- 
dicated by the slope of the free surface right at the well (cf. Fig. 3). 
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Fig. 5. Relation between fluid head and distance from outer sand 
face for various rates of flow. 


IV. THE QUANTITATIVE RESULTS 
Pressure distribution 


For purposes of application a quantitative formulation of the results of 
the investigation must be developed. From theoretical considerations it ap- 
pears that the phenomena should be described by an equation for the free 
liquid surface involving as parameters, the heights in the tank and well or 
the production rate from the system. Such a development is in fact outlined 
in §6. It will be found, however, that in attempting a rigorous analytical 
solution, it is impossible to carry through in detail a satisfactory treatment 
of the problem. Our procedure therefore has been that of analyzing the ex- 
perimental data to obtain an empirical formulation of the result. 
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As the most accurate data of the experiments were the readings on the 
manometer tubes, these were plotted to see if they implied any definite 
empirical relation. In Fig. 5 is plotted Z,°—Z,? against log. 7:/r, for fixed 
values of the production Q. Z, is the fluid height in the manometer tube 
farthest from the well at radius r; and Z,, is the manometer reading at radius 
r, taken at each of the points along the sector at which a manometer was 
attached. The production rate obtained for the pressure conditions repre- 
sented by each curve is shown. The curves plotted in Fig. 5 are calculated 
from the experimental data of Table I. It will be noted that within experi- 
mental error the curves are linear and pass through the origin. 
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Fig. 6. Relation between fluid head and distance from outer sand face for various rates 
of flow for determination of the constant A. 


In Fig. 6, the same data are plotted with ordinates Z,?—Z,” against the 
production rate Q for fixed values of log r;/r,. It will again be seen that the 
curves are linear and pass through the origin.* It follows at once that the fol- 
lowing relation must hold between (Z,?—Z,”), Q, and log r;/r, 


A(Z\> — Z,") = Q log. ri/r, (1) 


where: Z;, Z, are the fluid heads at the bottom of the sand at the points r;, 
r, respectively, and A is a constant. The value of the constant A for the 
twelve different sections of the sector represented by the various points r,, 
and at ten different production rates, is shown (Fig. 6) to average 0.078, 
with an average deviation from the mean within experimental error. This 
shows the validity of Eq. (1), and there remains the interpretation of the 
constant A in terms of known constants of the flow system. This determina- 
tion of A could not be made accurately in these tests because of uncertainty 


* The discrepancies occurring for large values of Q will be shown later to be due to the 
effect of the capillary layer. , 
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in obtaining the sand permeability and the unknown capillary effect and was 
therefore left for an additional experiment with the apparatus of Fig. 2 which 
will now be considered. 


Determination of the flow constant, A 


It will be seen from theoretical considerations which will be presented in 


§6 that in Eq. (1): 
A = rkyg. (2) 
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Fig. 7. Permeability determination, pure radial flow. Permeability 
k =264 cc/sec./cm*?/atmos./cm at 20°C. 


Considering the neglect of vertical velocities in the theoretical derivation 
we may reasonably expect, however, that the actually observed value of A 
will differ from rkyg by some appreciable amount which may be determined 
experimentally. 

Experimental data regarding the factor A obtained with the apparatus 
of Fig. 2 with the procedure which has already been described, will therefore 
be presented. 

Since the only unknown factor appearing in Eq. (2) is the permeability k 
of the sand, the first step is to determine its value. As mentioned in the de- 
scription of Fig. 2, this apparatus was suitable for observing flow rates under 
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pure radial flow conditions, that is, with the fluid levels at the input and out- 
put ends of the sector held above the sand height. 

Now the outflow from such a system is given rigorously (since the pres- 
sure has the role of a velocity potential) by the equation: 


Q = [2rkyeZ.(Z, — et l) log. r1/ Tw (3) 


where: Z,=sand height in cm; Z,;, Z~=fluid heights at input radius 7; and 
well radius 7, respectively. The value of & for the sand in place is therefore 
readily obtained. Fig. 7 shows the results of two different series of radial 
flows represented by the circles and dots respectively from which, we obtain 
the value of: | 

k = 264 cc/sec./cm?*/atmos./cm. 

A gravity flow was next observed holding the fluid level at the input end 
just below the sand height in order to eliminate as far as possible the capillary 
flow contribution. The pressure distribution and the rate of flow obtained 
are shown in Table III. 

TABLE III. Gravity flow. (Apparatus of Fig. 2.) 











Zi Z2 Z3 Zs Zs Zs aw Q 
(cm) cc/sec 
28.57 27.47 26.07 25.07 24.07 22.37 18.07 6.61 
r (cm) 39.0 30.0 21.0 16.5 12.0 7.5 3.0 


Sand height = 29 cm 








Applying the above data to the empirical relation previously obtained: 
A(Z,*? — Z,2) = Qlogeni/r, [cf. Eq. (1)] 


we obtain the following values for A, considering several values of r;/7r,,: 


loge 11/fw = 2.56 A, = 0.0346 
log. ri/re = 1.65 Ae = 0.0345 
loge7i/rs = 1.18 Az = 0.0329 
log. ri/ry = 0.862 Ag = 0.0305 


0.0331 average observed. 


It has been suggested that we might expect A to be given by rkyg 
modified by some additional factor. In fact, however, applying the value of 
the sand permeability k=264 cc/sec. determined experimentally (Fig. 7) 
and putting y=1.0, g=0.00097 atmosphere-cm/density and considering a 


15° sector: . 
A = (x X 264 XK 1.0 X 0.00097)/24 = 0.0335 


which is in closer agreement with the observed value than would be anti- 
cipated, considering experimental error. We are therefore justified in assum- 
ing that the factor rkyg is exact. 

The empirical gravity flow formula (1) may now be written 


QO = rkyg(Z,? — Z*)/log. ri/r (4) 

















TS 


mera 





FLOW OF LIQUIDS THROUGH POROUS MEDIA 103 


where: y=density of fluid, k=sand permeability (cc/sec./cm*/atmos./cm) 
Z,, Z are fluid heads in cm at r; and r respectively, g=acceleration of gravity 
=().00097 atmos.-cm per unit density, Q = fluid flow (cc/sec.). 

It is interesting to note that in form this formula is exactly that derived 
by neglecting vertical velocities. 


V. COMPOSITE RADIAL AND GRAVITY FLOW 
Analysis 


It is of practical interest to extend the discussion to include the case in 
which the fluid level at the input end of the system is greater than the sand 
height. Such a system corresponds to a gravity flow system in which the in- 
put fluid level is at the sand height with a superposed pressure at the outer 
boundary. Consideration of this composite flow system will suggest that 
analytically it may be treated as a purely radial flow and a gravity flow 
superposed. 
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Fig. 8. Diagrammatic representation of composite flow system. 


Fig. 8 shows diagrammatically the conditions assumed. Zo is the total fluid 
head at the outer boundary and Z, is the fluid head at the well. Z, is the 
sand thickness. The radial flow component is expressed as: 


O, = 2wkygZ.(Zo — Zs)/loge ri/rw- (5) 
The gravity flow component will be: 
Q, = whyg(Z.? — Zw)/l0ge 11/Tw. (6) 


The resultant composite flow will be the sum of the above components 


and is 
O., 5, = twkyg(2Z.Zo — Zs* ae Z~*)/ log. ri/ Vw (7) 


where: Z)=total fluid head in cm expressed in terms of the fluid flowing in 
the sand, Z,=sand height (cm), Z,, = fluid head at well (cm), 7; = outer radius 
(cm), 7»=well radius (cm). 

It is also of interest to obtain the equation showing the fluid head dis- 
tribution in the composite flow system. Again assuming a superposition of 








104 WYCKOFF, BOTSET AND MUSKAT 


the two types of flow, we readily obtain the expression 


(Zo — Z,) ‘ (Z,? sens Zw) — 
Pd = ee log. t/ Tw a Z;" _— ce — log. ri/r, (8) 


log. 71/Tw log. 1/T w 





where Z = hydrostatic head at radius 7. 


Experimental verification 


To test the pressure distribution indicated by the above analysis, Eq. 
(8), the apparatus of Fig. 1 having 20 manometers was selected as best suited 
for this purpose. The apparatus was modified by adding a tank extending 
above the sand level at the input end, thus permitting the water level at this 
end to be held above the sand level. In addition, after filling the sector with 
sand and water a layer of plaster of paris about 2 inches thick was poured 
over the sand top. This effectively sealed the top sand face for moderate 
pressure 

Table IV shows the comparison between calculated and observed fluid 
heads Z, and shows that, as far as pressure distribution is concerned, Eq. (8) 
represents the observed facts: 


yrs 


























TABLE IV. 
vient Z (calculated) | Z (observed) | _— Z (calculated) | Z (observed) 
(cm) (cm) (cm) (cm) 
6.4 0 | 0 | 55.8 39:9 40.0 
10.2 16.4 | 15.7 66.0 41.8 41.7 
15.3 23.5 | 22.8 | 76.2 43.0 43.1 ’ 
20.3 27.8 27.4 86.4 44.4 44.3 : 
25.4 30.6 | 30.6 | 96.5 45.6 45.4 
30.5 33.0 | 33.2 | 106.7 46.7 46.5 ; 
35.6 | 34.9 35.1 | 116.8 47.6 47.4 
40.6 36.5 | 36.8 | 127.0 48.3 48.5 
45.7 37.6 38.0 | 137.2 49.1 49.3 
50.8 390 | 38.9 || 147.3 49.7 50.0 
| 156.0 50.4 50.4 " 
Where the boundary conditions were: 
Zo = 50.4 cm, r, = 6.4 cm, 
Z,=40.5cm, re = 156.0 cm, 
Z, = 0.0cm, Q= 50 cc/sec. 


Verification of Eq. (7) expressing the flow rates is not obtainable from the 
experiment just cited because no independent value of the sand permeability 
k was obtainable. However, a new series of tests with the apparatus of Fig. 
2 with known sand permeability (Fig. 7) were made for this purpose. The 


essential observed data and the calculated Q’s are given in Table V. 
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TABLE V. Composite radial flow. 























Zo Le par vase Q Observed Q % error 
89.47 27.47 46.97 46.1 +1.9 
78.38 27.47 38.56 38.4 +0.4 
59.67 27.47 24.38 23.7 +2.9 
39 .37 27.47 8.99 8.6 +4.5 
70 .87 22 .57 36.07 36.1 —0.1 
47.57 22.59 18.41 18.5 —0.5 
78.37 17.17 44.56 , 44.2 +0.8 
58.47 By 29.48 29.0 +1.7 
44.97 17.17 19.25 19.3 —0.3 
74.37 12.57 43.32 43.5 —0.4 
59.57 12.57 32.09 31.5 +1.9 
48.17 12.57 23.46 23.9 —1.8 
68.17 6.97 40.05 40.8 —1.8 
50.87 6.97 26.94 27.0 —0.2 
35.55 6.97 18.32 16.4 —6.6 
61.47 0.17 35.61 36.1 —1.4 
56.47 0.17 31.82 Ce —0.9 
48.17 0.17 yee 260.7 —4.4 

Average error = —0.2% 








| 


Where: Z, = 29.0 cm, tkyg/24 = 0.0335, r,; =39 cm, 7, =3 cm. 

These data amply justify the conclusion that the composite flow formula, 
Eq. (7), gives accurate quantitative results and that considering either the 
pressure distribution within the system or the quantity output, this type of 
flow may be resolved into two discrete types; pure radial drive and gravity 
flow. In any practical application it is obvious that the superposed head im- 
posed at the outer boundary must be expressed as a fluid head in terms of the 
fluid flowing through the system. 


An approximate calculation of the capillary flow 


It has been mentioned that flow occurs in the capillary zone (see Fig. 4) 
and that the experiments from which the flow laws could be derived were 
necessarily conducted so that the contribution of the flow in this capillary 
layer was small compared with the total flow. Unless such precautions are 
taken the capillary flow will not be a negligible factor in the flow system. 
Thus in Fig. 3, flow No. 9, the capillary zone of 9.5 cm height is practically 
equal to the fluid head and as a result its contribution is by no means negligible, 
but is approximately 65 percent of the pure gravity flow. This is an extreme 
case but if calculations are made for other runs from the data given in Table 
II it will be found that there is a progressive increase in the discrepancy due 
to the increased proportion of capillary flow. It is therefore of some interest 
to obtain an approximately quantitative estimate of the flow occurring in this 
zone. In most practical cases this detail could be neglected but in many cases, 
and particularly in experimental work where the capillary rise cannot be 
made small without incurring serious deviations from pure viscous flow, an 
approximate calculation is of value. 

In Fig. 3 are plotted the fluid head distribution and the corresponding 
observed fluid levels in the sand for two sets of flow experiments. It should be 
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mentioned that the fluid heights in the sand cannot be obtained with any 
great accuracy even though a steady-state flow is maintained for several 
hours. However, these curves, which are typical of all of a large number of 
observations, show clearly the capillary layer and the discontinuity which ex- 
ists between the fluid levels in the well and in the sand face adjacent to the 
well. 

No attempt is made here to study the details concerning this discontinu- 
ity. Qualitatively, it is easy to see that vertical velocities must exist in the 
gravity flow system and that these vertical velocity components give rise 
to resistance forces in any vertical column which are proportional to the 
integrated vertical velocities in the column. Thus in Fig. 9 the actual fluid 
height Z, in the sand at any point, neglecting capillarity, is the observed fluid 
head Z, (measured at the bottom of the sand) +0Z,/kyg, where 3 is the aver- 
age vertical velocity of the fluid in the column. Unfortunately we do not know 
the velocity distribution in the system and hence cannot evaluate i Z,/kyg 
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Fig. 9. Diagrammatic representation of gravity flow system showing 
effect of capillary layer. 


which is the difference between fluid heads and fluid heights, existing every- 
where in the system. It is obvious, however, that as the fluid level in the well 
is lowered the vertical velocities increase and the discontinuity AZ, = (Z,,, 
—Z,-) between the fluid level in the well and in the adjacent sand face, for 
example, continuously increases. A rough analysis indicates that we might 
expect 4Z,, in the sand, to be at least approximately AZ/2 and the data of 
Fig. 3, together with other observations not shown, confirm this within the 
accuracy of experimental observations. Assuming that the fluid height in the 
sand is given by this relation we may proceed with the consideration of the 
capillary flow. 

Fig. 9 shows schematically an approximation to the actual flow system 
where Z; and Z,, are the observed fluid heads at 7; and +, respectively, Z; and 
Z.. are the fluid heights in the sand, and Z, is the thickness of the capillary 
zone. The surfaces indicated as P=0 are constant pressure surfaces. 

As already mentioned, we may expect the capillary zone to act like a 
radial flow syphon operating under a differential head 


(Z; po Zw) = (Z, = pm ; a 





a 
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Hence the radial flow in the capillary zone, neglecting end effects which are 
partially compensating, and assuming the capillary zone to be complete over 
the entire system, is 

QO. = [2rkygZ(Z, — Z~)/2|/loge ri/rw- (9) 


The purely gravity flow is: 
QO, = whyg(Z." — Zy")/log. ri/rw. (10) 
Hence the resultant flow is the sum of (9) and (10), or 
Qosre = whygl(Z: — Zw)(Zi + Zw + Z.)/log. ri/rw. (11) 


While this equation cannot be exact it should at least provide an ap- 
proximate correction for the capillary flow and through its use we will present 
additional experimental data on gravity flow which gives further support to 
the validity of the empirical Eqs. (4), (7), and (8). 

The data of Table VI cover observations on gravity flow using the ap- 
paratus of Fig. 2. The input and well fluid heads are given together with the 
associated production rates and the calculated values of production using the 
simple gravity flow formula, Eq. (10), and the calculated production cor- 
rected for a complete capillary zone (Eq. (11)), where: rkyg =0.0335 (meas- 
ured—Fig. 10). Z,=9.5 cm., log, 71/f%w = 2.56. 


























TABLE VI. 

Calc. Q Calc. Q Obs o7 

Run No. Z\ By Z, cc/sec. cc/sec. x /0 
Eq. (10) Eq. (11) ¢ erros 
1 31.47 18.07 29.0 8 .64* 8.61 + 0.4 
2 28.57 18.07 29.0 6.40 7.42 6.61 +16.6 
3 23.67 18.07 29.0 3.06 3.76 3.90 — 3.6 
4 15.37 0.17 29.0 3.09 4.98 4.95 + 0.6 
5 14.37 0.17 29.0 2.70 4.47 4.55 — 1.8 
6 9.57 0.17 29.0 1.18 2.38 2.50 — §.2 








* Calculated by Eq. (7)—composite flow, Z) =31.47 cm, and showing the value of & to be 
accurate. 

The capillary layer in run No. 2 was far from complete as will be noted, 
Z, being practically at sand height, Z,=29 cm. In this case the simple gravity 
formula gives a result only 3 percent low. The run is included here in order 
to show that the approximate correction may be made only in the case of a 
fairly complete capillary zone. No attempt has been made to handle the case 
of an incomplete zone. These data are intended to show the necessity of pre- 
caution in tests where appreciable capillary flow exists rather than to present 
a proof that Eq. (11) is accurate although it would appear that it is a fair 
approximation. 


VI. SomME THEORETICAL CONSIDERATIONS 


Strictly speaking the problem of gravity flow may be considered as a po- 
tential problem for the determination of the velocity potential 


& = k(p + gz), (12) 
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satisfying Laplace’s equation 
0°0/dx? + 0°/dy? + 0°/dz* = 0 (12a) 


and certain boundary conditions. These are: ®=constant at outside bound- 
ary, ®=constant at well surface. —06/dz=vertical velocity=0 at z=0, 
© =k(po+ygz); O&/dn=0, at the free surface whose normal is » and which 
is at the uniform pressure Pp. 

From this point of view the whole problem could be solved analytically 
by potential theory methods and the extended experimental investigation 
of the previous sections would be quite unnecessary. On the other hand, if 
we look at the boundary conditions, it seems that the problem is over-deter- 
mined since we specify both ® and its normal derivative at the free surface. 
And this discloses the fundamental difficulty of the problem which is that the 
free surface must be determined simultaneously with ® and in such a way that 
the last of the above boundary conditions are satisfied. The solution to this 
problem would evidently involve an integral equation of considerable com- 
plexity, although in two-dimensional cases the methods of conjugate functions 
could be applied, as has been done so elegantly by Hopf and Trefftz (refer- 
ence 3). 

To get a quantitative description of the problem of gravity flow it was 
therefore necessary to make the empirical investigation and analysis de- 
scribed in the previous sections. Nevertheless it is of interest to see what some 
simple assumptions as to the nature of the flow will give in the form of a 
qualitative description of the problem. We consider first the velocities at the 
free surface. 

Thus from the fact that ® is a velocity potential it follows that the 
velocity V, along the direction indicated by s which corresponds to the 
streamline h(s) is: 


V, = — k(Op/ds + ygdh/ds). (13) 


Since for the streamline defining the free surface, p is uniform, the velocity 


along the free surface is: 
V, = — kygdh/ds. (13a) 


The radial component is evidently : 
V. = — (kygdh/dr)/|1 + (dh/dr)?]. (13b) 


To get further we shall assume that the fluid moves towards the well in 

cylindrical shells so that the radial velocity in the interior of the sand is 

independent of the vertical coordinates and hence is equal to that at the free 

surface. It follows then that the total flux, Q, through the sand is given by: 
2arkyg(dh/dr) h(r) 


= — 2erV,h(r) = + 14 
e eee 1 + (dh/ dr)? (14) 





Taking Q as known, Eq. (14) leads to the following differential equation 
for the free surface h=h(r): 


(dh/dr)* — (lrkygrh/Q) dh/dr +1 = 0. (15) 





4 





od 
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Since this equation, too, being nonlinear, is difficult to solve explicitly we 
make the further assumption that the slope of the free surface, dh /dr, is small 
so that its square may be neglected. We then get: 


(2rkygrh/Q) dh/dr —1=0 (16) 
with the integral: 
h? — h,? = (O/rkyg) log r/ny. (17)° 


This equation is evidently identical with that established empirically in 
§4 except that here the / refers to the actual height of the free surface whereas 
empirically Eq. (17) will be correct if h refers to the fluid head as measured 
at the bottom of the sand. This relation, however, must be considered as a 
fortunate coincidence since a little inspection will show that the two assump- 
tions made in deriving Eq. (17) are definitely untenable and it is indeed sur- 
prising that the errors of these assumptions can be completely compensated 
by the simple change in the interpretation of h. 

Thus, since 

OV ,/dz = OV./dr 


our assumption of shell flow (0V,/dz=0) implies that the vertical velocity, 
V., is independent of r and hence must vanish since it vanishes at the external 
boundary. But without vertical velocities the free surface would be horizontal 
and the flow would have no gravity characteristics. Further, it is not difficult 
to show that not only will V, actually decrease as we approach the free sur- 
face, but even the total resultant velocity must decrease with increasing 2. 

That our neglect of (dh/dr)? in passing from Eqs. (15) to (16) is not strictly 
justified is still more obvious. For certainly near the well where the slope of 
the free surface approaches unity the first and last terms in Eq. (15) are of 
the same order of magnitude. The above derivation, which is equivalent to 
those in the literature except that it shows more clearly the nature of the as- 
sumptions made, can therefore be given no significance except as the present 
experimental study has established its form with regard to the pressure dis- 
tribution at the bottom of the sand. 

Nevertheless we can draw the following general conclusions from the 
above derivation. Thus, solving Eq. (13b), which is rigorously correct, for 


dh/dr we find 
dh kyg 1 gag i" 


& 
dr 2V, 2 v, — 


The plus sign before the radical would evidently correspond to a slope de- 
creasing with V, increasing. Since V, is negative (it decreases with increasing 
r) it follows at once that the maximum value of dh/dr is 1, or that the max- 
imum inclination that the free surface can have is 45°. When the inclination 
is 45° the radial velocity has its maximum 

| = kyg/2. (18a) 


| a 
r 
max. 


® This is essentially the derivation given by others. Handbuch der Physik. VII, 117, ex- 
cept that heretofore the assumption of shell flow and neglect of (dh/dr)? have always been 
made at the beginning. 
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It follows still more generally from Eq. (13a) that since the horizontal velocity 
component at the free surface of any gravity flow system is: 


Vi = — (yg sin 20)/2 


where @ is the inclination of the free surface to the horizontal, the maximum 
inclination of the free surface in any gravity flow system is 45° and the maxi- 
mum horizontal velocity component at the free surface is kyg/2, which is half 
of the maximum vertical free-fall velocity. This does not imply that the free 
surface inclination increases uniformly until a value of 45° is reached, since in 
the actual radial flow case due to the change in medium at the well boundary, 
the free surface cuts the well face horizontally giving an inflexion at some dis- 
tance from the well. Rather, the above results refer merely to the impos- 
sibility of obtaining inclinations exceeding 45° at any point of a gravity flow 
system. 

As to the slopes of the stream lines below the free surface the following re- 
marks may be made. Since the slope at any point is given by: 


m= V./V, 
we have: 
Om/dz = (1/V,)[0V./dz — m(OV,/dz) |. 

We have already seen that the radial velocity V,, decreases as the free 
surface is approached, so that 0V,/dz <0, considering the absolute values of 
the velocity components. Hence at such radii where the vertical velocities are 
increasing with the height the streamline slopes will increase as the free sur- 
face is approached and since the slope of the free surface cannot exceed 1, 
then neither will the streamlines have slopes exceeding 1. However, in the 
region near the well where the free surface undergoes an inflexion and its 
slope begins to decrease again the vertical velocity will also decrease so that 
for streamlines not far from the free surface it may well be that!® dV,/dz <0 
and finally that 0m/dz<0. Under such conditions it is possible that near the 
well the streamlines not much below the free surface may have slopes ex- 
ceeding 1. This does correspond to the observations described in §3 although 
it should be mentioned that the capillary layer would produce a similar effect 
and may also be effective in causing the sharp slopes of the streamlines as ob- 
served in the experimental study. 


VII. GENERALIZATIONS AND CONCLUSIONS 


The above experiments and discussion have referred to the problem of 
gravity flow in a radial system. Although in practical cases where a well of 
small radius is draining its surrounding sand, the assumption of radial flow 
must represent a close approximation to the actual flow system up to ap- 
preciable distances from the well center, it is of some importance to be able 
not only to test this assumption but even to be able to treat gravity flow prob- 
lems that definitely are not of the simple radial type. We shall therefore now 


10 This will in fact be necessary very near the well for small fluid heights in the well since 
then V)~kyg and since V,=0 at z=0 and small at the free surface it must be appreciably larger 
and even approach kyg as one recedes from the free surface and enters the sand body. 
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present a generalization of the above results which is felt to be the natural 
extension of our point of view as to the significance of the Eqs. (7) and (8) 
established empirically. 

Again as in the case of the radial flow problem, our generalization will 


‘ not be new in its analytical form. For as long ago as 1886 Ph. Forchheimer? 


presented the equation: 
0°2Z7/dx*? + 0°Z7/dy*? = 0 


as governing the general case of gravity flow problems. Here again, however, 
the Z’s are taken to refer to the fluid heights in the sand and the derivation of 
the equation is explicitly based on the assumptions of shell flow and the ne- 
glect of vertical velocities. These we have already seen to be physically un- 
justified as a region of fluid convergence is approached, and only fortuitously 
have they led in the case of radial flow to an analytical form empirically 
verified. Further, the computation of the flux at the well boundary would be 
totally in error if the above assumptions were logically applied since the 
actual discontinuity at the well face would give too large an area of outflow 
for the radial velocity computed by the above equation. 

We therefore propose to consider the problem from the following point of 
view. We observe only that the result of the empirical study on the radial 
flow system, which is: 





79 Z;° oe Z2* 4 79 
ZZ? = ictecta log ra + Z;? (19) 
satisfies the equation: 
0°Z*/dx*? + 0°Z7/dy*? = 0 (20) 


with the boundary conditions: 

Z=Zir=n; 

(21a) 

Z=Leir=fe 
where Z refers only to the fluid pressure head at the bottom of the sand. Since even 
a radial flow system corresponds to some degree of geometrical complexity 
we now generalize and postulate that Eq. (19) holds for every gravity flow 
system with a horizontal impermeable base with arbitrary boundary condi- 
tions. 

To get the total flux Q through a general gravity flow system we proceed 
as follows. Since at the high fluid level boundary the free surface slopes are 
necessarily small and the height of the free liquid surface actually coincides 
with the boundary value of the pressure head Z, we are justified in using the 
assumption of shell flow at this boundary and to identify the fluid height with 
the fluid head. We have then by D’Arcy’s law: 


pd fz ie ee ood (21 
= kyg Z—ds = — — ds 2 
‘ ” On 2 On 
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where :k, y, g have the same meaning as usual, ” is the external normal at the 
element ds of the boundary curve of the region, and the integral is extended 
about the whole of this boundary curve. Thus for the case of radial flow we 
have by Eq. (19): 


o mt f Zi 20 1 ag, rhe: — 2s") 
= a r 3 
r 


) 


2 log(ri/re) log(ri/re) 


which is identical with that given by Eq. (4) and verified empirically. 
For linear flow in a sand column of length L the solution to Eq. (20) 
would evidently be: 


Z? = (2,7 — Z2*)x/L+Z2 


so that 











Nm 
Nm 
— 


kyg Z,* — Z?* kyg Z,*> — Z2” 
gates ao eeese, ) 
L 2 L 


2 
where s is the width of the linear flow section. Eq. (22) has also been verified 
experimentally by measurements on a linear flow tank in a manner exactly 
the same as used for the radial flow experiments. This may be considered at 
least as a partial justification for the generalization of Eq. (20) for other than 
radial flow systems. 

If at the high fluid level boundary the liquid extends above the sand thick- 
ness we again have a composite flow case. The treatment then will follow ex- 
actly the lines of §5. Thus the gravity fluid heads Z, will be obtained by the 
solution of Eq. (21) with boundary conditions: 


Z, = Zz at boundary (2), 
Z, = Z, at boundary (1); 


where Z, is the sand thickness; then the nongravity component Z,, is to be 
found as usual from the equation 


8°Z,,/Ox? + d°Z,/dy? = 0 


i) 
w 
— 


with the boundary conditions 
Z, = Z, at boundary (2), 


; ; (23a) 
Z, = Z, at boundary (1), 


where Z, and Z» are the actual fluid heads at the boundaries (1) and (2) 
respectively it being assumed, of course, that 


Z,>2,>2Z2. 


The resultant fluid head distribution will then be: 


Z22=2Z,%+2,—-2.. 
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Similarly the resultant flow through the system will be given by: 


QO = vg [ (az, On),ds = kygZ, fz On)ds 


= (kyg/2) foz. On)ds + kygZ, fz. On)ds. 


These equations will suffice to determine everything of interest about any 
gravity flow problem in a homogeneous sand with a horizontal impermeable 
bottom boundary, with the exception of actual fluid heights in the sand 
throughout the system. In practical cases the effects of capillarity, which are 
not included in the above expressions for Q will be of negligible importance. 

The authors wish to acknowledge the assistance of Dr. H. T. Kennedy 
in some of the earlier experimental work on this problem. 
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An expression has been derived on the basis of the Maxwell-Boltzmann distribu- 

tion and Poisson's equation for the tangential surface conductance of an electrolyte 

in contact with a charged plane surface. One part of this expression is analogous to 

the result given by Smoluchowski for a Helmholtz double layer, while the second part 

arises from the movement of the ions in the diffuse ion cloud under the applied electric 

field. The theory shows qualitative agreement with the data of Briggs for the con- 
ductance of univalent chloride solutions at a cellulose surface. For all but the most 
dilute solutions, the Smoluchowski term is negligible and the surface conductance is 
proportional! to the charge density on the surface. The theory does not agree with the 
data for different anions and it is suggested that the discrepancy may be due to ab- 
sorption phenomena. There may be a movement of the adsorbed ions or of a diffuse 

ion cloud resulting from the presence of adsorption potentials at a distance from the 

interlace. 

oo the presence of a Helmholtz double layer at the surface of a 

solid in contact with an electrolyte, Smoluchowski, (1905), showed that 
there should be a surface conductivity due to the movement of the layer of 
charge in the electrolyte, when an electric field is applied parallel to the in- 
terface. Such conductivities have been reported by Stock (1912), McBain 
(1921, 1929), Fairbrother and Mastin (1924), Stamm (1926), Briggs (1928), 
and Gortner (1931, 1932), and his associates. It is now well known that the 
Helmholtz concept of a rigid double layer may be replaced by a diffuse ion 
cloud in which the concentration of ions in the neighborhood of a charged 
surface is given by the Maxwell-Boltzmann principle and the Coulomb 
forces, cf. Gouy (1910), Chapman (1913), Hertzfeld (1920), Gyvemant (1923), 
Stern (1924). 

Following Debye and Hiickel (1923), it will be assumed that we are deal- 
ing with a dilute solution of a strong electrolyte in which the average con- 
centration of the 7 ion is %; and its valence is 2;. Then, where the electric 
potential is y, the concentration of the 7 ion n; is given by the Maxewll- 
Boltzmann distibution function, 


n; = new tT = HAL — s:eh/kT), 1 


where ¢€ is the electronic charge, k is the Boltzmann constant and TJ is the 
absolute temperature. The free charge density is then 


p= Donze = — (Y/kT) dai. 
Substituting in the Poisson equation 
e 4rp 4rey ; 
div grad ¥ = — — = — ) fiz? = wy; 


D DkT 


1 The major portion of this work was done as an International Education Board Fellow 
at the Physikalisches Institut, Universitat Leipzig in 1928-29. 
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where 
x? = (4re/DkT) onze. 


For the semi-infinite one dimensional case with a plane at y=0, 
d*y/dy?= xy and yp = Ae + Ber, 
Since Y must remain bounded, B =0. When the surface charge density is ¢ 
(dy/dy)y-0 = — Ax = — 4n0/D 
and 


YW = (400/Dx)e-"". (2) 


By substituting in (1) 
o Anz An; 
a= £41 —- —-— cw) =€@ji- — ), 
xk DkT ni; 


If an electric field E is applied parallel to the plane, then because of the 
net charge density resulting from the unequal concentrations of the ions of 
opposite sign there is a body force on the liquid, 


where 


F = Ep = (— ED/4n) - (d*f/dy’). 
When the viscosity of the liquid is 7, for one dimensional laminar flow, the 
velocity u, is given by 
F + nd*u,/dy*? = 0 

and 
ED d*y du, 

— +9— =0. 
4nr dy? dy? 


Integrating, 
us = (ED/4xn)(W — Yo) 

since u,=0 when Y=y, and y=0. Then from Eq. (2) we have 
u, = — (Eo/nx)(1 — e~*”) 


which is the motion of the liquid due to the body force. 
Since this body force is the result of the viscous drag of the ions on the 
liquid we must also consider the velocity of the ions relative to the medium, 


u; = Eez;/p; 


where p; is the friction coefficient of the 7 ion, so the total velocity of the ion 
relative to the plane is 
u;’ = Us + Uj. 
The current density 


An; 
J= DYoufnze= Yo(u, + ui) (1 - —) HiZie 


nj 
= Do uit ie ~ > uAn zie 
+ Do uinzie - douAnize. 
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Subtracting the current density due to the bulk conductance X of the elec- 
trolyte we have Lf;zie=0 and LA,z;eu;=X and the surface conductance cur- 
rent density is 


: Eo Eo 4Are® nz? 
AJ =J— Ey = —-on(l — oo) - —:-—— or ) 
nk k DkT pi 


AX = —dy= ie | - e (3) 
The first term is the counterpart of Smoluchowski’s result and is due to the 
net charge carried by the motion of the solvent while the second term, due 
to the motion of the ions through the solvent, apparently has not been con- 
sidered except by Mooney (1931). 
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Fig. 1. AA, the surface conductivity in reciprocal megohms vs. o, the surface charge density in 
e.s.u./cm?’, from data of Briggs on cellulose. 





Since either o or Yo must be known, the only data available with which to 
test Eq. (3) are for the surface conductance of packed cellulose fibers—of 
which those of Briggs (1928) are representative. These data for the univalent 
chlorides have been corrected for the ionic strength of the distilled water by 
Abramson (1932), who showed that there is a general parallel between 
o*/x and AX as predicted by Smoluchowski. It is found, however, that this 
term becomes smaller with increasing salt concentration since x increases 
faster than o? and that for a uni-univalent salt in which the mobilities of the 
two ions are quite different, the second or mobility term is considerably more 
important. On the assumption that for all except the lowest concentrations 
the Smoluchowski term would be negligible, ¢ has been plotted against AX 
in Fig. 1. The linear relationships at the higher concentrations are quite 
striking in their qualitative agreement with the theory. It is unfortunate that 
it is not possible to make an accurate check on the correctness of the magni- 
tude of the specific surface conductance as given by the theory because of the 
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complex nature of the surfaces in such a diaphragm of packed cellulose fibers. 
An approximate calculation based on an assumed cell length of 1 cm gave 
the sum of the perimeters of the transverse sections for all fibers of such a 
diaphragm as about 104 cm. This does not seem an impossible value. Assum- 
ing a reasonable value for Yo, Mooney (1931) calculated AX for glass to be 
0.02 of the conductance given by McBain (1929). Since White, Urban, and 
Van Atta (1932) extrapolate a value 0.03 of that of McBain’'s, they agree 
rather well with Mooney. 

There is, however, a very marked disagreement between the theory and 
experiment. Since the mobility term is proportional to the differences of the 
mobilities of the two ions, it should not contribute to the surface conductance 
of KCl. This is most decidedly not the case,—in fact the effect of Cl- seems 
too large for all of the salts. The approximation involved in the expansion of 
the exponential in Eq. (1) is not serious enough to explain the discrepancy 
when the maximum value of Yo is 13 mv. The interionic forces in the electro- 
lyte and their modification in an electric field have not been considered. These 
effects should be treated by the methods of Onsager (1927) and Debye and 
Falkenhagen (1928), but Wagner (1924) showed the former effect to be 
negligible for the negative adsorption at an air-electrolyte interface. The as- 
sumption of a semi-infinite electrolyte at a plane interface may be replaced 
by one of a tube of electrolyte bounded by a cylindrical interface. If the radius 
vo of this tube is small so that nyo < <1, the approximations to the Bessel’s 
functions of the solution lead to the result that Yo, o and AX are independent 
of x. This is not-in agreement with the data. 

By comparison of the surface conductances of the various electrolytes, 
at a given surface charge density, with the ionic mobilities, it is found that 
the effect of the Cl- ion is fairly constant and considerably more important 
than is indicated by Eq. (3). This suggests that either the adsorbed charge or 
the adsorption forces may play a role in these phenomena. If the adsorption 
forces are not effective at distances greater than molecular or ionic dimen- 
sions, then there is the possibility of a net motion of the adsorbed ions under 
the influence of a tangential electric field as well as a two dimensional Brown- 
ian movement, cf. McBain and Salmon (1921), Volmer (1922). In this case, 
the surface conductance due to the adsorbed ions would be 


where ; and p; are the “concentration” and “mobility” of the ion on the sur- 
face. This hypothesis will explain quite well the experimental observations 
for uni-univalent electrolytes, where the form of the ¢ vs. concentration curve 
suggests a preponderate adsorption of the negative ion, cf. Abramson and 
Miiller (1932). It is also in agreement with the observation of Briggs (1928) 
that AA¥0 when y¥o=0 for certain polyvalent electrolytes. If ions of both 
signs are adsorbed such that for some concentration ¢ =0, both the Smolu- 
chowski and the mobility conductances will vanish, while this adsorption 
conductance is not necessarily zero. Langmuir (1918), Gyemant (1923), 
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Stern (1924) and Hiickel (1928) have assumed that the sphere of action of 
adsorption forces is of the order of molecular dimensions, but according to 
Trillat (1926) it can be shown that for liquids at an unsaturated interface 
these forces may be effective at a considerably greater distance. The cor- 
responding adsorption potential would cause the formation of a diffuse ion 
cloud near the interface in much the same manner as the electric potential 
considered above. This adsorption ion cloud would then give rise to an added 
surface conductance without the assumption of mobility of ions adsorbed at 
the interface. The analysis of this case by the Maxwell-Boltzmann distribu- 
tion would require for adsorption potentials the counterpart of Poisson's 
equation in electrostatics—since the inverse square law of force action would 
not be expected to be immediately applicable. 

In conclusion it might be pointed out that such adsorption factors as 
these considered may play an important part in other electrokinetic phe- 
nomena and that careful experiments of the type performed by Briggs but 
of a more absolute character are greatly to be desired. 

I wish to express my thanks to Professor P. Debye and to Dr. H. A. 
Abramson for their interest and assistance in this problem. 
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Interaction of Light and Sound Waves 


N A recent paper by Debye and Sears at- 

tention is directed! to a very pretty experi- 
ment on light scattering in liquids. Originally 
suggested by L. Brillouin® in 1922 from the- 
oretica! considerations, it has recently been 
carried out in connection with experiments 
made on the Raman effect in liquids. 

It is a familiar fact that in a homogeneous 
medium the propagation of light takes place 
very much as it does in free space; i.e., in 
“straight lines” and with uniform velocity. 
(That the velocity of propagation is less than 
in free space is not of importance from this 
particular point of view.) On passing the 
boundary separating two such individually 
homogeneous regions the light is subject to a 
change, however, and generally speaking is 
partly transmitted (refracted) and partly re- 
flected, so that the discontinuous variation in 
the optical properties at the boundary sepa- 
rating the two media leads to a “breaking up” 
of the light, to which the general term scatter- 
ing can be applied. 

Naturally, the transition boundary between 
two media is not really discontinuous, but can 
be thought of as a thin layer within which the 
properties of the one medium vary rapidly to 
those of the other. If we wish to treat the 
boundary as though it were perfectly sharp 
we can do so by considering it as the ideally 
limiting case of a very thin but finite transi- 
tion layer. 

Similarly, a medium whose properties vary 
in a continuous manner from point to point 
can be treated as the limiting case of a con- 
glomerate of small homogeneous regions with 
discontinuous boundaries. As a physical illus- 
tration we recall the familiar fact that sand 
of which the individual particles are very 
minute has many of the properties of a liquid, 
for example, that of flowing. If one grinds the 


! P. Debyeand F. W. Sears, Proc. Nat. Acad. 
Sci. 18, 409 (1932). See also E. H. L. Mayer 
and W. Ramm, Phys. Zeits. 33, 270 (1932). 


?L. Brillouin, Ann. de Physique 17, 88 
(1922). 


sand grains finer and finer into a powder the 
sand behaves more and more like a liquid. 
What we wish to see from this for our present 
purpose is that in traversing such a conglom- 
erate of minute particles (such as a “sand” 
made of glass particles) light would be re- 
fracted and reflected at each one of the bound- 
aries between particles, so that it would be 
gradually dispersed or “scattered.” Even in 
the limiting case in which the sizes of the par- 
ticles become infinitesimal this scattering ef- 
fect persists, even though the medium can be 
finally treated as being continuous (but of 
course not uniform). 

The upshot of the matter is that inhomo- 
geneities in a substance traversed by light, 
will lead to a scattering of the light. It is im- 
material how the inhomogeneities are pro- 
duced, whether by mechanical stresses, non- 
uniform heating or in any other way. We may 
recall in passing the extensive use which is 
being made of this fact in the testing of 
stresses in bridges and the like by constructing 
glass models through which light is passed. 
The inhomogeneities produced by the strains 
in the glass produce visible colorations which 
can be readily studied. 

The basic idea of the experiments treated 
by Debye and Sears is that if an otherwise 
optically transparent medium is transmitting 
sound waves then the fluctuations of pressure 
within the sound wave will produce inhomo- 
geneities in the density of the transmitting 
medium of just the type which we have been 
discussing, so that if light is passed through 
the medium it should be scattered. We are 
thus justified in speaking of an interaction of 
light and sound, or of a scattering of light by 
sound. These terms hardly do the phenome- 
non justice but at least they give succinct 
expression to it. 

This idea of Brillouin’s opens up a number 
of suggestive lines of thought. For example, if 
we consider the simplest case of a plane sound 
wave of wave-length /, then the distance be- 
tween a region of maximum pressure and one 
of minimum pressure in the wave is //2, and 
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since these regions are separated in space by 
regular intervals, they present an array simi- 
lar to a grating. Can one then expect to find 
an action on the light similar to the action of 
a grating? Will one find interference phenom- 
ena associated with fluctuations in intensity 
of the light “reflected” from this grating? 

The answer given by theory is in the affirm- 
ative, and this has been verified by the experi- 
ments. In the arrangement of Debye and 
Sears the “sound waves” were generated arti- 
ficially in a liquid by immersing a piezoelectric 
crystal in it. On setting the crystal into oscil- 
lation its motion was transmitted to the liq- 
uid. The liquid was contained in a vessel with 
transparent walls through which light could 
easily be transmitted, and the scattered light 
examined as it emerged from the other side. 
The light scattered at different angles with 
the direction of the primary beam was ex- 
amined in order to look for the interference 
patterns of which we have spoken, and which 
are due to the grating character of the sound 
waves. These patterns were readily found. 

It is of particular interest to note that all 
of the scattered light does not have the same 
wave-length as the incident light. This is due 
to the fact that the sound waves in the liquid 
are travelling and not standing waves. If we 
use the simplified idea of the scattering proc- 
ess as being simply a reflection of the light 
from a set of moving parallel planes (or mir- 
rors) then we must expect that there will be 
a Doppler effect giving rise to scattered light 
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of frequency »[{1+2(v/v’) sin (@/2)] where vo 
is the frequency of the incident light, v’ the 
velocity of the light in the liquid, v the veloc- 
ity of the sound waves in the liquid and @ the 
angle between the directions of the incident 
and scattered light. Since there is always 
some scattering due to dust particles, motes, 
and the like in the liquid, (and also some scat- 
tering due to the Raman effect of individual 
molecules) without change of frequency of the 
light, one must expect to find light of three 
frequencies in the scattered radiation, and of 
course if one could measure the change in fre- 
quency of the scattered radiation accurately 
a determination of v, the velocity of sound in 
the liquid, could be made. This has been veri- 
fied experimentally and preliminary measure- 
ments of this nature are reported by Debye 
and Sears. 

Actually the experimental results are some- 
what more complicated and require a more 
detailed theory than we have been able to 
give here, principally due to the finite size of 
the actual apparatus which precluded the 
possibility of setting up purely harmonic 
plane waves in the liquid. The authors state 
that the scattered light can easily be made of 
sufficient intensity for projection on a screen, 
which of course cannot be done with the feeble 
scattered radiation in the Raman effect, and 
which would serve to give an audience a vivid 
picture of the process of light scattering. 


E. L. Hitt 


Quantum Mechanics and Optics 


HROUGHOUT the whole period of de- 
velopment of quantum mechanics con- 
siderable interest has attached to the analogy 
with mechanics, an 


analogy which can be roughly stated in-the 


optics and Newtonian 


form of a proportion: 


Quantum mechanics: Newtonian mechanics: : 
Physical optics: Geometrical optics. 


The proportion is of course imperfect, but is 
indicative of the formal analogy. The relation- 
ship between geometrical optics and Newton- 
ian mechanics was discovered years ago by 
the eminent geometrician and logician, Sir 
William Rowan Hamilton, but the comple- 
mentary relationship between quantum me- 
chanics and physical optics was not discovered 


until about 1924 when L. de Broglie developed 


his ideas on the wave nature of matter. The 
latest link in this chain of reasoning is fur- 
nished in a paper by W. Pauli! who has an- 
alyzed the relationship existing between the 
relativistic wave equation of Dirac and the 
dynamics of the special theory of relativity. 
During the early period of development of 
quantum mechanics there was a great deal of 
discussion concerning the role played by the 
Planck constant h. It was observed to be a 
general rule that whenever a formula found 
by means of quantum theory did not contain 
h, then that formula could be (and in most 
cases already had been) developed directly 
from classical electron theory: and if, on the 
other hand, h appeared explicitly in the for- 
mula then the classical theories gave either 
bad answers or else no answers at all to the 
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problem. This made it seem pretty certain 
that there must be a close relation between 
the equations of the “true” quantum theory 
and those of the classical theory, associated 
in some manner with the limit h—0. So far 
as the nonrelativistic wave equation of Schré- 
dinger is concerned, the method for making 
this comparison was developed independently 
by G. Wentzel and L. Brillouin. 

Before delving into this problem let us ex- 
amine some of the physical elements of the 
problem. In the Newtonian dynamics, say for 
a single particle, one proposes to obtain equa- 
tions stating the coordinates of the particle as 
a function of time; i.e., x(t), y(t), 2(t). The 
method for doing this is to solve the differen- 
tial equations for the motion, assuming that 
the analytical expressions for the components 
of the force are known; to complete the solu- 
tion one must specify the way in which the 
particle is started off, that is, its initial posi- 
tion and velocity. When all of these things 
have been stated the solution of the equations 
of motion will give one the path or trajectory 
followed by the particle, and they will give in 
addition all information concerning the way 
in which the particle traverses this path— its 
velocity at any point, its time of arrival at any 
point, etc. In other words, we get a complete 
time table for the touring particle. This beau- 
tifully simple idea of a point travelling along a 
definite curve in a perfectly definite manner 
can easily be extended to more complex sys- 
tems by considering the motion of a single 
“representative” point in a (mathematical) 
space having as many dimensions as the sys- 
tem has degrees of freedom. 

But instead of thinking of a single particle 
being started off from a given point at some 
definite instant we can think of a very large 
number of particles all having the same total 
energy but being started off from a variety of 
points. The particles are all supposed to move 
in the same “external” field. of force without 
influencing each other. We shall then have the 
picture of a swarm of “touring” particles all 
travelling along their individual paths. Since 
the potential energy of a particle may be con- 
sidered as a definite function of position, the 
speed of the particles depends only on the 
particles’ positions (and of course on its total 
energy which is the same for all particles by 
hypothesis). 

The analogy with geometrical optics can 
now be made by the following observation, the 
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proof of which is given by the Hamilton- 
Jacobi dynamical theory. Suppose we set up 
in the field of force a surface, the form of 
which is arbitrary (but which for simplicity 
we may visualize as a portion of a sphere). 
Let us consider a swarm of particles all start- 
ing at the same instant from different points 
on the surface with the velocities appropriate 
to those points and travelling away from the 
surface in the same direction. The paths of 
these particles will then form a family of 
curves orthogonal to the given surface. If we 
now imagine the family of surfaces orthogonal 
to this family of curves (and of which the orig- 
inal surface is one member) drawn in, then the 
proof referred to shows that all of the par- 
ticles move together in such a way that at any 
instant they all lie together on a single mem- 
ber of this family of surfaces. Or, if we 
imagine a surface drawn at every instant in 
such a way that the set of particles all lie in 
it, then this surface, as it travels along with 
the particles, continuously adjusts itself so 
that it “fits” some member of the family of 
surfaces. The analogy with geometrical optics 
is now clear—the family of surfaces are the 
analogs of the surfaces of constant phase for 
the light, while the paths of the particles; i.e., 
the orthogonal trajectories to the surfaces are 
the analogs of the rays of the light. If at this 
point one takes the trouble to adjust the equa- 
tions of motion of the particles, in particular 
the force equation, then it is possible to de- 
velop a quite sensible theory of photons to 
replace the mass particles with which we com- 
menced. In this way geometrical optics, com- 
bined with dynamical theory leads at once to 
a corpuscular theory of light. 

This is all very well, but we know that ge- 
ometrical optics fails to give any description 
of such important phenomena as scattering, 
diffraction, dispersion, etc. Just so does New- 
tonian dynamics fail to give an account of the 
analogous phenomena for matter waves as we 
now understand them; and just as there is no 
bridge for the gap between Fermat's theorem 
which lies at the basis of geometrical optics 
and the optical wave equation which forms 
the foundation of physical optics, so there is 
no bridge for the gap between Newton's equa- 
tions of motion and Schrédinger’s equation. 
There is the analogy, and fortunately this has 
been sufficient to give us a start. 

We know that the differences between the 
predictions of physical and geometrical optics 
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become ironed out in the case of light with 
small wave-length. Practically, of course, this 
takes us into the regions of x- and y-rays when 
neither theory really applies, but formally the 
two sets of equations merge into each other if 
we take the limit \—>0, where A is the wave- 
length of the light. If we look for an analogous 
degeneration of quantum mechanics into 
Newtonian mechanics we find it in precisely 
the same way when we take as the wave- 
length of a mass particle the de Broglie ex- 
pression \=h/mv, m being the mass of the 
particle, v its velocity, and h Planck’s con- 
stant. 

But now we see the meaning which is to be 
attached to the formal limit 4-0 which we 
spoke of at the beginning of this note, and 
which seems a little awkward since h is a per- 
fectly definite and finite quantity! The answer 
is that the proper “physical” limit is not h-0, 
but A—0 which can be approached by con- 
sidering m— © or better v-@. For high-ve- 
locity particles then, those distinguishing 
features which separate the quantum me- 
chanical description of their behavior from 
that of Newton’s dynamics disappear. 

The only fly in the ointment now is that it 
is known that for high-velocity particles one 
should replace the Newtonian dynamics again 
by the dynamics appertaining to at least the 
special theory of relativity. And as for the 
wave equation of Schrédinger, we guess that 
Dirac’s relativistic equation is the proper first 
substitute. 

The partial justification that this procedure 
is in the right direction constitutes one of the 
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most interesting points in Pauli’s paper,! in 
that he actually carries through the demon- 
stration that the dynamical equations of rela- 
tivity theory are the zero’* order approxima- 
tion to Dirac’s equation. The demonstration 
proceeds along the lines of the previous treat- 
ment of Wentzel and Brillouin which we have 
already mentioned, and which in essence is 
an expansion of the solution of the quantum- 
theory equations as a series in powers of / 
(although not in a simple Taylor’s or power 
series). 

In the first approximation involving h, 
which is that in which deviations from the 
“geometrical optical” properties of the rela- 
tivistic dynamics become apparent, Pauli 
finds that the peculiar effects due to the spin 
of the electron also appear. By studying the 
behavior of the electron in a one-dimensional 
field of force he has been able also to study 
more carefully the so-called Klein Paradox 
according to which a mass particle can change 
the sign of its mass, an electron thus behaving 
effectively like a proton of mass 10?’ g, for 
the detailed consideration of which we must 
refer to the original article. 

It is especially interesting to notice in these 
investigations how a small variation in the 
formal appearance of a set of equations can 
change so radically their properties and the 
associated physical interpretation. 


E. L. Hitt 


1W. Pauli, Helvetica Physica Acta 5, 179 
(1932). 








